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We introduce and analyze a novel metallic phase of two-dimensional (2d) electrons, the Roton 
Fermi Liquid (RFL), which, in contrast to the Landau Fermi liquid, supports both gapless fermionic 
and bosoiiic quasiparticle excitations. The RFL is accessed using a re- formulation of 2d electrons 
consisting of fermionic quasiparticles and hc/2e vortices interacting with a mutual long-ranged sta- 
tistical interaction. In the presence of a strong vortex-antivortex (i.e. roton) hopping term, the RFL 
phase emerges as an exotic yet eminently tractable new quantum ground state. The RFL phase 
exhibits a "Bose surface" of gapless roton excitations describing transverse current fluctuations, has 
off-diagonal quasi-long-ranged order (ODQLRO) at zero temperature {T — 0), but is not supercon- 
ducting, having zero superfluid density and no Meissner effect. The electrical resistance vanishes as 
T with a power of temperature (and frequency), R{T) ~ T'' (with 7 > 1), independent of the 
impurity concentration. The RFL phase also has a full Fermi surface of quasiparticle excitations 
just as in a Landau Fermi liquid. Electrons can, however, scatter anomalously from rotonic "cur- 
rent fluctuations" and "superconducting fluctuations" . Current fluctuations induced by the gapless 
rotons scatter anomalously only at "hot spots" on the Fermi surface (with tangents parallel to the 
crystalline axes), while superconducting fluctuations give rise to an anomalous lifetime over the 
entire Fermi surface except near the (incipient) nodal points ("cold spots"). Fermionic quasiparti- 
cles dominate the Hall electrical transport. We also flnd three dominant instabilities of the RFL 
phase: an instability to a conventional Fermi liquid phase driven by vortex condensation, a BCS- 
type instability towards fermion pairing and a (non-pairing) superconducting instability. Precisely 
at the instability into the Fermi liquid state, the exponent 7 saturates the bound, 7 = 1, so that 
R{T) ~ T. Upon entering the superconducting state the rotons are gapped out, and the anomalous 
quasiparticle scattering is strongly suppressed. We discuss how the RFL phase might underlie the 
strange metallic state of the cuprates near optimal doping, and outline a phenomenological picture 
to accomodate the underdoped pseudo-gap regime and the overdoped Landau Fermi liquid phase. 



I. INTRODUCTION 

Despite the appeal of spin-charge separation as an un- 
derpinning to superconductivity in the cuprates, there 
are seemingly fatal obstacles with this approach. Ever 
since Anderson's initial suggestion— of a spinon Fermi 
surface in the normal state at optimal doping, there have 
been nagging questions about the charge sector of the 
theory. The concept of the "holon", a charge e spin- 
less boson, was introduced in the context of the doped 
spin liquid state^, and was presumed to be responsible 
for the electrical conduction. In the simplest theoretical 
scenario, spin-charge separation occurs on electronic en- 
ergy scales, thereby liberating the electron's charge from 
its Fermi statistics. This idea has been actively investi- 
gated over the past 15 years - see Ref. IH for a review. 
A pervasive challenge to this perspective, however, is the 
difficulty of avoiding holon condensation and supercon- 
ductivity at inappropriately high temperatures. In addi- 
tion, recent theoretical work, which has elucidated the 
phenomenology of putative spin-charge separated states, 
has led to further conflicts with observations. One class 
of theories have shown how spin-charge separation can 
emerge from a superconducting phase by pairing and con- 
densing vortices.'* Following this work, a Z2 gauge theory 
formulation greatly clarified the nature of fractionaliza- 
tion of electronic (and other) quantum numbers.^ It has 
become clear that a necessary requirement for true spin 



charge separation in two dimensions is the existence of a 
"vison" excitation with a gapALSiS The vison is perhaps 
most simplest thought of within the vortex pairing pic- 
ture, as the remnant of an unbound single vortex. If these 
ideas were to apply to the cuprates, one would expect this 
gap to be of order the pseudogap scale ksT* . Unfortu- 
nately for spin-charge separation advocates, experiments 
designed to detect the vison and measure its gapifi*ii have 
determined an unnaturally low upper bound of approxi- 
mately 150if for the vison gap in underdoped YBCO. Is 
this the death knell for spin-charge separation? 

In a very recent paper focusing on the effects of ring- 
exchange in simple models of bosons hopping on a 2d 
square lattice^^, we have identified a novel zero tem- 
perature normal fluid phase - (re-)named the "Exciton 
Bose Liquid" (EBL). In the EBL phase boson-antiboson 
pairs (ie an exciton) are mobile, being carried by a set 
of gapless collective excitations, while single bosons can- 
not propagate. The resulting quantum state is "almost 
an insulator" , with the d.c. conductivity vanishing as a 
power of temperature, cr(T) ^ with a > 1. This is in 
contrast to the "strange metallic phase" in the optimally 
doped cuprates, which is "almost superconducting" 1^ 
with an extrapolated zero resistance at T = 0, as if it 
were a superconductor with Tc = 0. This phenomenol- 
ogy suggests the need for a non-superconducting quan- 
tum phase in which the vortices are strongly immobilized 
at low temperatures. 
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Motivated by this, we revisit the Z2 vortex-spinon field 
theory of interacting electrons^, in which the hc/2e vor- 
tices and the spinons have a long-ranged statistical in- 
teraction mediated by two Z2 gauge fields. Rather than 
gapping out single vortices while condensing pairs (which 
leads to a spin-charge separated insulator)^, we would 
like to find a quantum phase in which the vortices are 
gapless but nevertheless immobile. To this end, we add 
an additional "vortex-ring" term to the earlier vortex 
field theory. This term is efiectively a kinetic energy for 
vortex-antivortex pairs, that is for rotons. To access the 
limit of strong roton hopping requires a further reformu- 
lation of the Z2 vortex-spinon field theory, replacing the 
Z2 gauge fields by U{1) gauge fields. Similar U{1) vortex- 
fermion field theories have been explored in Refs.[l3- The 
resulting U{1) vortex-spinon formulation is tractable in 
the limit of a very strong roton hopping, and describes 
the "Roton Fermi Liquid" (RFL) phase, a novel metallic 
ground state qualitatively different than Landau's Fermi 
liquid phase. 

Both the Z2 vortex-spinon field theory developed in 
Refs. 4 5] and our U{\) vortex-fermion field theory, de- 
scribed in Section II below, are constructed in terms of 
operators which create the excitations of a conventional 
BCS superconductor: the Bogoliubov quasiparticles, the 
hc/2e vortices and the collective plasmon mode. This 
choice of "basis" , however, does not presume that the sys- 
tem is necessarily superconducting at low temperatures, 
and indeed we intend to employ such a formulation to de- 
scribe non-superconducting states. When superconduc- 
tivity is present at low temperatures, the formulation will 
also be employed to describe the "normal" state above 
Tc- In these approaches, the Bogoliubov quasiparticle is 
electrically neutralized^^, and the resulting "spinon" exci- 
tation transported around an hc/2e vortex within an or- 
dinary 2d BCS superconductor, acquire a Berry's phase 
of TT. Within the Z2 and U{1) vortex-spinon formulations, 
one introduces spinon creation operators, fl^, where we 
let r denote the sites of a 2d square lattice and with a 
the spin. Vortex creation operators are also introduced, 
conveniently represented in a "rotor" representation as 
e'^, which live on the plaquettes of the 2d lattice. The 
vortices are minimally coupled to a gauge field, living on 
the links of the dual lattice. The "flux" in the gauge 
field describes charge density fluctuations on the original 
lattice sites, and for example encapsulates the plasmon 
mode inside the superconducting phaseii. Finally, the 
long-ranged statistical interaction between the spinons 
and vortices is incorporated by introducing two Chern- 
Simons Z2 or U{V) gauge flelds. The important new el- 
ement in the present paper is the inclusion of a roton 
hopping term. As we shall see, the RFL phase is readily 
accessed within the U{V) formulation when the magni- 
tude of this roton hopping term is taken signifigantly 
larger than the single vortex hopping strength. 

Here, we briefly summarize the main results estab- 
lished in the following Sections. In Sec. II the Z2 vortex- 
spinon field theory formulation of Ref. (sl is recast in 



terms of a lattice Hamiltonian. Via a sequence of exact 
unitary transformations on the Hamiltonian, we demon- 
strate that it is possible to exchange the Z2 Chern-Simons 
gauge fields for their C/(l) counterparts. Within a La- 
grangian representation of the resulting U{\) Hamilto- 
nian which we employ throughout the paper, it is possible 
to choose a gauge for one of the C/(l) Chern-Simons fields 
so that the "spinon" is recharged, and has finite overlap 
with a bare electron. (In Apd.IdIwc show that this gauge 
choice can effectively be made at the Hamiltonian level, 
and construct a Hamiltonian theory in terms of vortices 
and fermionic operators which carry the electron charge 
and have a finite overlap with a bare electron.) 

Initially, in Sec. HI, we ignore the fermions entirely, 
and focus on the bosonic charge (or vortex) sector of the 
theory. A "spin-wave" expansion valid in the presence 
of a large roton hopping term, leads to a simple theory 
which is quadratic - except for a single vortex hopping 
term. Dropping this vortex hopping term then leads to 
a soluble harmonic theory of the "Roton-Liquid" (RL) 
phase. In addition to the gapless 2d plasmon, the RL 
phase is shown to support a "Bose surface" of gapless ro- 
ton excitations. We compute the Cooper pair propagator 
in the RL phase, and show that it exhibits off-diagonal 
quasi-long-ranged order (ODQLRO) at zero tempera- 
ture, but not a Meissner effect. The RL phase exhibits 
a high degree of "emergent" symmetry - the number of 
vortices on every row and column of the 2d dual lattice is 
asymptotically conserved at low energies. This symmetry 
implies that the harmonic "fixed point" theory of the RL 
phase has an infinite conductivity at any temperature. 

In Sec. IV we study the legitimacy of the approxi- 
mations used to arrive at the harmonic RL theory, fo- 
cussing first on the neglected vortex hopping term. We 
show that for a range of parameters the vortex hop- 
ping term is "irrelevant" , scaling to zero at low ener- 
gies whenever its associated scaling dimension satisfies 
At, > 2. Nevertheless, at finite temperatures vortex hop- 
ping leads to dissipation, giving a resistance which van- 
ishing as a power law in temperature, R{T) ~ with 
7 — 2A„ — 3 > 1. A "plaquette duality" transformationi^ 
allows us to next address the legitimacy of the initial 
"spin-wave" expansion, used to obtain the harmonic RL 
theory. Of paramount importance is the presence of a 
term in the dual theory which hops a "charged" quasi- 
particle excitation, a term not present in the harmonic 
fixed-point theory. We find that the "charge" hopping 
process is irrelevant over a range of parameters - approx- 
imately the complement of the range where vortex hop- 
ping was irrelevant - implying stability of the RL phase. 
When relevant, on the other hand, the "charge" quasi- 
particle condenses, leading to a superconducting ground 
state. 

The fermions are re-introduced back into the theory in 
Section V, where we argue for the stability of the Bose 
surface of rotons and the 2d plasmon in the presence 
of a gapless Fermi sea of fermionic quasiparticles. We 
denote the corresponding phase by the Roton Fermi Liq- 
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uid (RFL). The gaplessness of the fermions is somewhat 
surprising, and deserves some comment. Indeed, it is 
in sharp contrast to the gapped nature of the quasiparti- 
cles in both the superconducting phase and Zi fractional- 
ized insulator, in which the spinous experience a BCS-hke 
"pair field" . Tlic cause of this difference is the existence 
of gapless single vortex excitations (and fluctuations) in 
the RFL, which according to our analysis leads to the 
"irrelevance" of the fermion pairing term. Crudely, be- 
cause the bosonic pair field exhibits only ODQuasi-LRO 
rather than ODLRO, there is no average pair field felt by 
the quasiparticles, and hence no gap. In this sense, the 
RFL is in fact closer to a Fermi-liquid state than it is a 
superconductor. 

Section VI is devoted to an analysis of the proper- 
ties of the RFL phase. We study both the longitudi- 
nal and Hall conductivities, and find that the dissipa- 
tive electrical resistance vanishes with a power of tem- 
perature, i?(r) ^ T'^ with 7 > 1, similar to the be- 
havior without fermions present. But the fermions arc 
found to dominate the Hall response, leading within a 
naive Drude treatment to an inverse Hall angle varying 
as G^'^ ^ 1/{tjT^), with Ty the fcrmionic quasiparti- 
cle transport lifetime. Due to the presence of the "Bose 
surface" of gapless rotons, electrons at finite energy w 
experience anomalous scattering, not present in a Lan- 
dau Fermi liquid. Specifically, quasiparticles scatter due 
to rotonic "current fluctuations" and "superconducting 
fluctuations" , which contribute additively to the electron 
decay rate. The former gives rise to especially strong 
electron scattering at "hot spots" - points on the Fermi 
surface with tangents parallel to the axes of the square 
lattice. At such hot spots the associated electron decay 
rate varies with an anomalous power of energy, cj'^^+^^Z^, 
for 1 < 7 < 2. The decay rate from superconducting fluc- 
tuations is present everywhere along the Fermi surface 
except near "cold spots" at the incipient d-wave nodal 
points. This contribution grows strongly with decreas- 
ing energy/temperature, although it has a smaller over- 
all amplitude than the current fluctuation contribution. 
Upon entering a superconducting state, the rotons - gap- 
less within the RFL phase become gapped out, and all 
anomalous scattering is strongly suppressed. 

Finally, in Section VII we briefly discuss possible con- 
nections of the present work to the cuprates. We suggest 
that the RFL phase might underlie the unusual behavior 
observed near optimal doping in the cuprates, in par- 
ticular the "strange metal" normal state above Tc. A 
scenario is outlined which also incorporates the pseudo- 
gap regime and the conventional Fermi liquid behavior 
in the strongly overdoped limit. 



II. THE MODEL 

We are interested in electrons hopping on a 2d square 

lattice, with electron creation operators cj^. Here, the 
sites are denoted (in bold roman characters) as, r = 



xix + X2y, where xi,X2 arc integers, Xi = x, X2 — y 
are unit vectors along the x and y axes, and a =t, i de- 
notes the two spin polarizations of the electron (such a 
spin index will be distinguished from Pauli matrices 
by the lack of any superscript). 

A. Z2 Chargon-spinon formulation 

We begin by formulating the electron problem in spin- 
charge separated variables using the Z2 gauge theory 
Hamiltonian. We emphasize that this formulation does 
not imply that spin-charge separated excitations are de- 
confined, and indeed this formulation correctly describes 
the low-energy physics of conventional confined phases as 
well. 

The Z2 gauge theory is most readily formulated in 
terms of a charge e singlet bosonic chargon br,bl, a 
neutral spin- 1/2 fermionic spinon ,fra,.f}cn ^^'^ Ising 
gauge field (Paidi matrix) <yj{r) residing on the link be- 
tween sites r and r + It is convenient to use a rotor 
representation for the chargons, br = e""^"", 6^6^ = '^r, 

with [(/'rj'^r'] = "iSr.r'- 

The Z2 Hamiltonian is conveniently expressed in terms 
of the Hamiltonian density, = which 

in turn is decomposed into a bosonic charge sector, a 
fermionic sector and a gauge field contribution, Ti.z2 — 

Ti-c = -^cXlo"|(r) COs{dj(j)r - Aj{r)) + UcjUr - pof, (1) 

j 

j □ (r+w) j 

j 

-te J2 /rV*,ae*^^^'^'-^^-^"'»/r. + h.C. (3) 
3 

Here, po = 1 ~ is the electron (charge) density with x 
measuring deviations from half-filling. Throughout the 
paper, dj with j = 1,2 denotes a discrete (forward) spa- 
tial lattice derivatives in the xi and X2 directions, for 
example, di4>r — dx4'r — f^r+x — i'r- We have included 
an external (physical) vector potential Aj (r) in order to 
calculate electromagnetic response and to include applied 
fields. The Hamiltonian Ti^ describes the dynamics of the 
chargons hopping with strength, tc^ which arc minimally 
coupled to the Z2 gauge field. The dynamics of the gauge 
fields is primarily determined from Hg, the first two terms 
of which constitute the standard pure Z2 gauge theory 
Hamiltonian. The "magnetic" contribution involves the 
plaquette product, 

n a^ = af(r)af(r + y)a|(r)a|(r + x), (4) 

□ (r-l-w) 

which is the Z2 analog of the lattice curl. Here we have 
defined, w = (l/2)(i; -|- y), and r -|- w denotes the cen- 
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ter of the plaquette. We have also included an additional 
contribution bi- linear in , which in the dual vortex rep- 
resentation below will become a "roton" hopping term. 
In the fermion Hamiltonian, H ^ ^ , we have defined the an- 
tisymmetric matrix e^a' = i'^^aui i E^nd take Aj = (— 1)^ A, 
which describes a nearest-neighbor pair field with d-wave 
symmetry. Apart from the first two terms familiar to afi- 
cionados of the Z2 gauge theory^, we have included a less 
exotic bare electron hopping amplitude ig- We will pri- 
marily be interested in the limit that the spinon hopping 
strength is significantly larger than the electron hopping 
strength, tg ^ te- 

In most of the analysis of this paper, we will consider 
the limit of small spinon pairing Aj 0. This can be 
justified either by the assumption Aj <^ ts, or by the ir- 
relevance in the renormalization group (RG) sense, which 
will occur occur in some regimes. If strictly A^ =0, both 
fermion (spinon) number and boson charge are conserved, 
and in principle may be separately fixed. However, for 
Aj 0, even infinitesimal, this is not the case. Instead 
the spinous will equilibrate in some time that diverges 
as Aj —^ but is otherwise finite, and the system will 
choose a unique fermion density to minimize its (free) 
energy. We will return to this point in the U{1) formu- 
lation in Sec, men 

The full Hamiltonian above has a set of local Z2 gauge 
symmetries, commuting with each of the local operators, 

Cl = (-1)"^+"^ H a^ (5) 

+ (r) 

where n( = f^^f^^ is the fermion density and the Z2 
lattice divergence is defined as: 

l[a^^l[a^{r)a^{r-±,). (6) 

+ (r) 3 

Physical states are required to be gauge invariant, which 
is specified by the set of local constraints: = 1. This 
is the Z2 analog of Gauss law for conventional electro- 
magnetism. 

The connection between the Z2 gauge theory and a 
theory of interacting electrons, is most apparent in the 
limit that t^, is taken to be much larger than all other 
couplings. In this limit the electron creation operator 
is equivalent to the product of the chargon and spinon 
creation operators, cl^ — h].fl^. Indeed, when ty ^ 00 
the Z2 "electric" field becomes frozen, ctJ « 1, and the 
gauge constraints then imply that on each site of the lat- 
tice the sum of the chargon and spinon numbers, n^ + nl, 
is even. Moreover, for large ty, the chargon and spinon 
hopping terms are strongly suppressed, and can be con- 
sidered perturbatively. Upon integrating out the gauge 
fields, one will thereby generate an additional electron 
kinetic energy term with amplitude of order tctg/ty. A 
brief discussion is given in Appendix IXI 

In what follows, we will study the Z2 gauge theory 
more generally, away from the large ty limit. Of interest 



is the electron Greens function, 

^e(ri,Ti;r2,T2) = -(r^Cri^(ri)c]:^^(T2)). (7) 
We will express the electron operators as, 

Cro- = fcr/ro-, (8) 

which is exact as t„ — > cx), but more generally should 
be sufficient to extract the universal low energy and long 
length scale behavior of the electron Green's function. 
We will also be interested in correlation functions involv- 
ing the Cooper pair creation and destruction operators, 

St, Br, with Br = {b^Y = 6-2'-/"'. 

B. Z2 Vortex-spinon formulation 

In what follows, it will prove particularly convenient 
to work with vortex degrees of freedom, rather than the 
chargon fields. To arrive at such a description, we use the 
U(l) duality transformation^^, in which the dual vari- 
ables sit naturally on the 2d dual lattice. We denote 
the sites of the 2d dual lattice by sans serif characters 
as r = xix. + a;2y + w, with w = (l/2)(x -I- y) and inte- 
ger xi,X2- The duality transformation itself defines two 
conjugate gauge fields [ai(r), ej(r')] — iSijSry, where 

rir = -ey9iao(r - w), (9) 

TT 

di4>r = 7reyej(r + w — Xj). (10) 

Here, as defined, due to the discreteness of the vari- 
ables, aj{r) takes on values that are integer multiples of 
TT, while ej{r) is a periodic variable with period 2. This 
transformation is faithful provided the constraint 

(V.e)(r)^^9,e,(r-x,) = (mod2), (11) 

j 

or equivalently 

^ e'^(^-^i') = 1 (12) 

is imposed at every site r of the dual lattice. Rewriting 
the charge Hamiltonian, one has 

V-c = -tc^ cr,f (r) cos(7reyej(r + w - Xj) - yli(r)) 

i 

+^{€ijdiaj - wpoY. (13) 

Conventional hc/2e superconducting vortices are com- 
posites of a "vison" (topological excitation in ct|) and a 
half-vortex in (p. To describe them, we perform a uni- 
tary transformation to a new Hamiltonian with new 
constraints C° , 

Hz2-^U^Hz,U, r^UTU, (14) 
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with the unitary operator, 

U = exp[^ eijaijr + w - x,)(g|(r) - 1) 



nK(r))"^(a2^(r))^^, (15) 



with w = (1/2) (x — y) = w — y. The transformed con- 
straints are, 



The fields e^^' and TV,, have the interpretation of vor- 
tex creation and number operators, as can be seen by 
tracking the circulation defined from the original char- 
gon phase variable One finds canonical commutation 
relations, 



+ (r) 

Cf = (-l)(v-'^(0 Yl a'- = 



(16) 
(17) 



□ (r) 



Under this unitary transformation, 

Uc = -<c^cos(7rej(r) + ejfeAfc(r + w)) 
j 



(18) 



=-t.^a^-(r)cos(a,(r))-if(-l)(^"^W +Hftl9) 

j 

Here we have defined, CTi(r) — erf (r -t- w),CT|(r) = (r — 
W), and have used Eq. (|17ll . The transformed "roton" 
hopping term becomes, 

= -K,a?(r)af(r-t-X2)cos(5ya,(r))-f(x ^ y). (20) 

The fermion Hamiltonian is almost unchanged in the dual 
vortex-spinon representation, 

i 

-t,Ea|(r)e'(-^-W-^^M)/;^-^.,/r.-f/..c., (21) 

j 

the changes appearing only in the electron hopping (the 
last term). Here we have defined. 



ei(r) =eyej(r + w-xj). 



(22) 



Notice that ifej couples "like a gauge field" to the spinous 
in the final electron term. 

To arrive at the final Z2 vortex-spinon theory, we 
split the electric and magnetic fields into longitudinal 
and transverse parts, aj = a' -I- a*, e, = -I- e*, with 



- , ^j, 

V • a* = V • e* = 0, eijd,a\ = eijdie] = 0. For fu- 
ture purposes, we note that the longitudinal part of Cj 
is related to the transverse part of Cj and vice versa, i.e. 
e'^*(r) = eije*^'(r-|-w — fj). It is convenient then to solve 
the constraint for the longitudinal fields, 

4(r) = -d,e,, (23) 
(V-e*)(r) = N,. (24) 



(25) 



At this stage, A^r is a periodic variable with period 2, 
and —dj9r + aj{r) is constrained to be an integer multiple 
of TT. It is convenient to soften the latter constraint, and 
in order to respect the uncertainty relation implied by 
Eq. (|25|) . at the same time relax the periodicity of N^. 
Formally, this is accomplished by replacing, 

-tc'^cos{nej+ejkAk) Const. +^ y^(ej+ejfcAfc/7r)^, 
J j 

(26) 

with Ui, ~ tcTT^, and adding a term to the Hamiltonian, 



2„y^cos(2aj6'- 2a*). 



j 



(27) 



The constraint is recovered for large t2v , but we will con- 
sider a renormalized theory in which t2v may be consid- 
ered a small perturbation. 

It is convenient to regroup the longitudinal contribu- 
tion to He along with the terms in Hg and into vor- 
tex "potential" and "kinetic" terms, Hm and i?^^- We 
thereby arrive at the final form for the Z2 vortex-spinon 
Hamiltonian, 



Hz-, — H, 



pi 



H^ 



N 



H 



kin 



H 



Z2 
f ' 



(28) 



with the fermion Hamiltonian li.f'^ given in Eq. (|21|l and 
with. 



^[4(r) + e,,4(r + w)]^ 



H, 



N 



n 



kin 



Z2 



Z2 



(29) 

)V{r-r'), (30) 



(31) 



Here B^ = eijdiAj{r — w) is the physical flux through 
the plaquette of the original lattice located at the site r 
of the dual lattice. In the "plasmon" Hamiltonian, Hpi, 
an implicit sum over j = 1 , 2 is understood and we have 
defined a (bare) plasmon velocity, vq, as Vq = 2uctc = 
2uci„/7r^. Inside the superconducting phase, this Hamil- 
tonian describes the Goldstone mode - or sound mode 
- and can be readily diagonalized to give the disper- 
sion, ujpii^) — I'd Z]j[2sin(fcj72)]^ with \kj\ < tt in the 
first Brillouin zone. Since the electron charge density is 
given by, ^eijOiUj, one can readily include long-ranged 
Coulomb interactions which will modify the plasmon at 
small k. 
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In Hf^ above we have set K = Q, dropping hence- 
forth the term 6Hn{K) = - K Y,^{-1)'^' , smce it wih 
not play an important role in the phases of interest. 
The vortex- vortex interaction energy V{r) is the Fourier 
transformation of the discrete inverse Laplacian opera- 
tor, V^-i(k) EE/C2(k),with/C2(k) =^^.2(l-cosfcj),and 
has the expected logarithmic behavior at large distances, 
V{r) ^ ^ln(|r|). The vortex kinetic energy, H^^„, is a 
sum of three contributions - a single vortex hopping term, 

= (^) cos(5,^, - 4), (32) 

a pair- vortex hopping term 7i^^ given in Eq. (|27|) , and a 
"roton" hopping term, 

= ~Kral{r)al{r+SLi)cos{/^.^ye,-d^al{r))+{x ^ y), 

(33) 

where we have defined, 

A,y0,^Ay,e,= J2 (-l)'^+''^r+eix+e.y. (34) 

61,62=0,1 

Notice that H^'^ hops two vortices, originally at sites of 
the dual lattice on opposite corners of an elementary 
square, to the other two sites. Equivalently, this term 
can be interpreted as hopping a vortex-antivortex pair 
on neighboring sites (ie. a vortex "dipole" or more sim- 
ply a roton) one lattice spacing in a direction perpendic- 
ular to the dipole. Such a roton is a 2d analog of a 3d 
vortex ring, and in a Galilean invariant superfluid (such 
as 4 — He) vortex rings propagate in precisely this man- 
ner. Henceforth we shall refer to this process as a "roton 
hopping" process. 

The above Z2 vortex-spinon Hamiltonian must be sup- 
plemented by the two gauge constraints, which from 
Eqs. (|15I17|I and (|24|l can be cast into an appealingly 
simple and symmetrical form: 

ci = (-1)"^ n = 1' (35) 

□ (r) 

Cf = {-!)''' II a^ = I. (36) 

□ (r) 

These constraints correspond to an attachment of a Z2 
flux in cr^ and to the vortex number parity and the 
spinon number parity, respectively. Since the spinous are 
minimally coupled to cr^ and the vortices to , this im- 
plies a sign change upon hopping a spinon around a vor- 
tex - or vice versa. Indeed, if the partition function for 
the vortex-spinon Hamiltonian (without Ti.^^ and with 
te — 0) is expressed as an imaginary time path inte- 
gral with the Z2 constraints in Eqs. 1^ imposed, 
the resulting Euclidean action becomes identical to Eqs 
(109)-(113) in Ref. 5]. 

It is instructive to obtain explicit expressions for the 
electron and Cooper pair creation operators in terms of 



the dualized vortex degrees of freedom. From Eq. ® 
we can directly obtain an expression for the Cooper pair 
destruction operator, = e^^^'f''^ , as 

00 

r 

where the Y[ symbol here denotes a product along a semi- 
infinite "directed" string running on the links of the orig- 
inal lattice, originating at r and terminating at spatial 
infinity, with dr' the unit vector from the point r' along 
the string to the next point. In terms of ej (rather than 
Cj), the product contains a factor of exp(±27riej(r)) for 
each link of the dual lattice that crosses the string, tak- 
ing the positive/negative sign for directed links crossing 
the string from right/left to left /right proceeding from r 
to 00. We will use the above notation when possible to 
present precise analytic expressions for such strings. The 
path independence of the string is assured by the second 
gauge constraint, Cf — 1. Since the unitary transforma- 
tion, U in Eq. (|15() , commutes with e^'^"^^ ^ this is the 
correct expression for the Cooper pair operator within 
the Z2 vortex-spinon theory. 

An expression for the electron operator, Cra = b^fra, in 
the dual vortex-spinon theory can be extracted by more- 
over re-expressing the "chargon" operator = e^^'^" as 
a string, 

00 

6r = n ^""'''"^ = S.,ort{r)S4r). (38) 

r 

For later convenience, we have here decomposed this ex- 
pression into a piece depending on the vortex configura- 
tions through the longitudinal "electric" field, e^, and a 
contribution depending on the smooth part of the phase, 
(/), through the transverse field, e*: 

00 00 
5™rt(r)-n^"'^'""^'; S4r)^Y[e'^~'>'^. (39) 

r r 

But unlike the Cooper pair operator, the "chargon" op- 
erator transforms non-trivially under the unitary trans- 
formation in Eq. (|15|) : 

00 

r 

now including a factor of cr| (r) for each link of the string. 
Again, the path independence is guaranteed by the gauge 
constraint, = 1. The final expression for the electron 
operator within the dualized Z2 vortex-spinon field the- 
ory follows simply as, 

Cra ^W'CraU ^bjra- (41) 

As discussed at length in Ref. , the Z2 vortex-spinon 
formulation is particularly well suited for accessing spin- 
charge separated insulating states. Specifically, when 
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pairs of vortices hop around they "see" an average gauge 
flux of, 2tijdiaj — 2TrpQ. Thus at half- filling with pq — 1, 
vortex-pairs are effectively moving in zero flux and at 
large pair-hopping strength, t2v — *■ oo, will readily con- 
dense - driving the system into an insulating state with 
a charge gap. In the simplifying limit with vanishing 
single vortex hopping strength, ty = 0, all of the vor- 
tices are paired, (—1)^ = 1, and the Z2 gauge constraint 
in Eq. (|^ reduces to Hn ''"'^ — 1- The "vison" excita- 
tions (plaquettes with JJ^ — —1) are gapped out of 
the ground state, and the spinous, being minimally cou- 
pled to cr^, can propagate as deconfined excitations. The 
charge sector supports gapped but deconfined chargon 
excitations, which can be viewed as topological defects 
in the pair-vortex condensate. 

But as we shall see below, to access the new Ro- 
ton Fermi Liquid phase requires taking the strength of 
the roton hopping strength large, and the Z2 formula- 
tion proves inadequate. To remedy this, we introduce 
in subsection C below, a new U{1) formulation of the 
vortex-spinon field theory. As we shall demonstrate, 
the U{1) and Z2 vortex-spinon formulations are formally 
equivalent, and by a sequence of unitary transforma- 
tions it is possible to pass from one representation to the 
other. Care should be taken when considering operators 
which transform non-trivially under the unitary opera- 
tions related different representations, however. A third 
dual vortex formulation involving electron (rather than 
spinon) operators is briefly discussed below in Appendix 
IdI The Hamiltonian in this "vortex-electron" formula- 
tion is equivalent under a sequence of unitary transfor- 
mations to both the Z2 and U{1) vortex-spinon Hamil- 
tonians. 

To establish these equivalences, it is convenient to 
"choose a gauge" in the Z2 theory. As detailed in Ap- 
pendix^ it is possible to unitarily transform to a basis 
in which the Z2 gauge fields are completely "slaved" to 
the vortex and spinon operators, and can be eliminated 
completely from the theory. Specifically, in the chosen 
gauge the x— components of both and are set to 
unity on every hnk of the lattice: af{r) = af{r) = 1. As 
we shall see in subsection C below, the U{1) gauge fields 
in the U{1) vortex-spinon formulation can be similarly 
enslaved. Remarkably one arrives at the identical "en- 
slaved" Hamiltonian in both cases, thereby establishing 
the formal equivalence between the two formulations. 

C. U(l) Vortex-spinon formulation 

In the U{1) formulation of the vortex-spinon field the- 
ory, the Pauli matrices , which live as Z2 gauge fields 
on the links of the original and dual lattice, respectively, 
are effectively replaced by exponentials of two U{1) gauge 
fields: o^{r) exp[iaj{r)] and cr|(r) exp[if]j{r)]. 
These two U{1) gauge fields are canonically conjugate 
variables taken to satisfy, 

K(r-x,),/3,(r')] =i^e»,<5'(r-r'), (42) 



with r = r -I- w. For two "crossing" links these com- 
mutation relations imply that the exponentials, e*",e''', 
anticommute with one another, [e'°',e'^]- = 0. 

1. U{1) Vortex-spinon Hamiltonian 

The full Hamiltonian for the U{1) vortex-spinon field 
theory takes the same form as the Z2 vortex-spinon 
Hamiltonian in Eq. f^ . 

H = Hpi+HN+ Hkrn +Hf, (43) 

with Hpi and given as before in Eq. (f^ . Only the 
vortex kinetic terms and the fermion Hamiltonian are 
modified. Once again the vortex kinetic energy terms 
are decomposed into single vortex, pair-vortex and roton 
hopping processes: 

'Hkin ='Hv + 'H2V + 'Hr. (44) 

Since the vortices in the U{1) formulation are minimally 
coupled to the U{\) gauge field, aj(r), each of these three 
terms will be modified from their Z2 forms. Specifically, 
in terms of the associated Hamiltonian densities we have, 

T-ty = -ty ^ cos(9j 61 - a* + a^), (45) 
'H2y = -t2y cos(2c)j6l - 2a) + 2a,), (46) 

Ur - -y cos[A,y0 - d,{al - ay)] + {x ^ y), (47) 

with l^xyQr defined in Eq. (|^ . 

The Hamiltonian density for the fermions in the U{\) 
formulation is given by 

3 

fra' + h.c] (48) 

In the U{1) formulation, the (average) density of spinous, 
{fa-fa) = i'^^ij'^i'^j) (^^ foUows from Eq. H54|l ') is taken 
to be equal to the (average) charge density, {■^(-ijdiaj) . 
A new element, not present in the Z2 fermion Hamilto- 
nian, Tif^ in Eq. H21(l . is the "string operator", Sr- The 
string operator is given as a product of e^^ running along 
directed links of the original lattice from the site r to 
spatial infinity: 

00 

5r = PJe"^^-''"-^-. (49) 

r 

As we shall discuss below, once we restrict the Hilbert 
space to gauge invariant states other choices for the 
"path" of the string are formally equivalent. 
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In addition to global spin and charge conservation, 
the full U{1) vortex-spinon Hamiltonian, Ti. above, has 
a number of local gauge symmetries. To fully define 
the model we need to specify the set of gauge invari- 
ant states that are allowed. Associated with each of the 
"Chern-Simons" fields, a and (3, is a U{1) gauge symme- 
try. Specifically, the full Hamiltonian is invariant under 
independent gauge transformations: 



^ Oij{r)+djXr, (50) 

and 

fro- ^ fva^ ^7 

/3,(r) ^ /5,(r)-fa,Ar, (51) 

for arbitrary real functions, Ar and Xr, living on the orig- 
inal and dual lattices, respectively. The corresponding 
operators which transform the fields in this way are, 

g^(Xr) = e-''Sr>^'[^-i^-«'ft(^— )1, (52) 

and 

g^(A,) = ^^E.K[fLfr.-iu,d,c,Ar-^)]_ (53) 

Both of these operators commute with the full Hamilto- 
nian H. The U{1) sector is specified by simply setting 
Gv ~ Gf ^ ^ ioT arbitrary Xr and A^. From Eqs. H52I53|) 
this is equivalent to attaching tt flux in the statistical 
gauge fields a and (3 to the spinons and vortices: 

eijd.iaj{r-w) = nf^^fra] e,j(9i/3j(r-w) = nN^. (54) 

Notice that this is simply the U{1) analog of the Z2 flux 
attachment in Eqs. (|35I36|I . and implies the same sign 
change when a spinon is transported around a vortex 
or vice versa. The only difference is that in the U{1) 
formulation the phase of tt is accumulated gradually when 
the spinon is taken around the vortex, whereas the sign 
change in the Z2 theory can occur when the spinon hops 
across a single link. The formal equivalence of the Z2 
and U{1) formulations will be established below. 

As we shall see, the "smearing" of the accumulated tt 
phase change, makes the theory in the U{1) formulation 
eminently more tractable. The one notable complication 
is the square of the "string operator", which in the U{1) 
sector is a non-trivial function of e'^^^ along the string, 
rather than equaling unity as in the Z2 sector. How- 
ever, it is worth emphasizing that within the U{1) sector 
of the theory, the value of the operator Or = is in- 
dependent of the chosen path. Specifically, consider two 
(unitary) string operators, denoted Oi, ©2, with different 
paths running from the same site r off to spatial infin- 
ity. The "difference" between the two string operators, 
Oi^02, is a product of e^''^ around closed loops. But due 
to the U{1) gauge constraint in Eq. eijdi(3j = nN, 

this product is an exponential of the total vorticity Ntot 



inside the closed loops, Oi^02 = exp(2niNtot)- Since 
the vorticity is integer, one deduces that the string oper- 
ator is indeed path independent, Oi — ©2. 

It will prove useful to obtain expressions for the elec- 
tron and Cooper pair creation operators within the C/(l) 
vortex-spinon formulation. The Cooper pair operator has 
the same form as in the Z2 vortex-spinon formulation, 
given explicitly in Eq. 1)37(1 . but the electron operator is 
modified in a non-trivial way. As we shall check explicitly 
below, the electron operator within the U{1) formulation 
involves a string depending on both the dual "electric 
field" , ej , as well as the statistical gauge field, (3j : 

00 

r 

with ei(r) defined in Eq. (|22|l . The path independence 
follows from the condition in Eq. ((24|l . together with the 
second U{1) gauge constraint in Eq. ifSH) above; 

TTdjCjir - Xj) = eijdi/3j{r ~ w). (56) 

This implies an equality between the longitudinal "elec- 
tric field" and the transverse statistical field: 

/3Kr)=7re!,(r-fw); /3* (r) = -e^r - w), (57) 

or 

A* = -^el (58) 

and enables the electron destruction operator to be re- 
expressed as, 

00 

Cra = fr.S^{r)l[e'^'^'"''^ (59) 

r 

with S^{r) defined in Eq. (|39() . Similarly, the string oper- 
ator that enters into the U{1) fermion Hamiltonian, Hf 
in Eq. ((48|l . can be written as, 

Sr = S,ort{r)Y[e'^'^'"''^ . (60) 

r 

Upon combining the above two equations, and recalling 
the identity for the "chargon" operator in the original 
Z2 theory, br = 5„ort (i")'?^ (r) in Eq. lp!H|l . imphes that, 
Srfra = blcra ■ Consequently, "spinon pairing" terms in 
the Z2 gauge theory are seen to be equivalent to the usual 
Bogoliubov-deGennes form: 

[Srf eaa' fra fra' = Ecrcr' Cro-Cro-' , (61) 

with BI the Cooper pair creation operator. For d-wave 
pairing, the pair field lives on links, and a similar identity 
obtains: 

i/3j(r)rc 12 /• f ,f , — r ^~ e ,r ,lei2\ 
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Here we have introduced a bond-Cooper pair operator, 
^L+x, = ^J^U' (63) 

OO 

bl = Y[e-'''^''^< . (64) 

r 

Note that in the un-transformed chargon variables of the 
Z2 gauge theory formulation, -Br.r+x^ = ''■|('^)^i^i+Xj ■ 

Recalling the discussion in Sec. Ill Al it is still the case 
that for Aj = 0, when this pairing term is absent, the 
fermion number is conserved, and naively may be cho- 
sen arbitrarily. However, in the limit Aj 0, which we 
consider here, this conservation is weakly violated, and 
only the total charge is conserved. The physics at work 
is clear from Eq. H61|l : for non- vanishing Aj, quasiparti- 
cle pairs and boson pairs are interchanged, and the two 
charged fluids come to equilibrium. Thus in what follows, 
we should choose to divide the charge density amongst 
the fermions and bosons in such a way as to minimize the 
total (free) energy. This division will therefore shift as 
parameters of the Hamiltonian are changed. How it does 
so is crucial to the ultimate low energy physical proper- 
ties of the system, as is clear from Eqs. (03- (03, which 
show that the vortices experience an effective "flux" pro- 
portional to the difference of the total charge density 
{eijdiaj/n) and the fermionic density Uf = eijdiUj /t:. As 
the fermion density is varied, the effective flux seen by the 
vortices changes. Significantly, in the limit t^, — > 00, vor- 
tex hopping dominates the energetics, and is minimized 
when the fermion density equals the total charge density, 
Uf = tijdiajj-K. This naturally recovers the Fermi liq- 
uid phase (App. ^ by binding charge e firmly to each 
fermion, fully accommodating all the electrical charge. 

In the rest of the paper we work exclusively within 
the U(X) vortex-spinon formulation, which is particularly 
suitable for extracting the properties of the Roton Fermi 
liquid. Before embarking on that, we first establish the 
formal equivalence between the two formulations by en- 
slaving the C/(l) gauge fields. As detailed in Appendix 
O it is possible to unitarily transform to a gauge with 
V • a = V • /^^''"'^ = 0. In this gauge, both aj and /3j are 
"enslaved" , being fully expressible in terms of the spinon 
and vortex densities, nl and iV,., respectively. Moreover, 
the enslaved ?7(1) Hamiltonian is found to be identical to 
the enslaved Zi Hamiltonian obtained in Appendix ^ 
and the enslaved expressions for the electron operators 
also coincide. 

Having thereby established the equivalence between 
the Z2 and [/(I) vortex-spinon field theories, in the 
remainder of the paper we choose to work exclusively 
within the C/(l) formulation, employing the Hamiltonian 
H defined in Eq. (|43|) . together with the gauge con- 
straints in Eq. (|54() . In practice, it is far simpler to work 
within a Lagrangian formulation, where the gauge con- 
straints can be imposed explicitly within a path integral, 
as detailed in the next subsection. 



2. Lagrangian for U{1) Vortex-spinon theory 

In order to impose the C/(l) gauge constraints in 
Eq. 154|) on the Hilbert space of the full vortex-spinon 
Hamiltonian, H, we will pass to a Euclidean path inte- 
gral representation of the partition function. The associ- 
ated Euclidean Lagrangian is readily obtained as a sum 
of three contributions, 

S = J dT[H + Lb + L,on], (65) 

with Lb involving the generalized coordinates and con- 
jugate momenta. 

Lb = ^ [iN.doO, ~ -(3i{r)€ijdoaj{r + Xi)] 

r— r+w 

+ H{r)doa]{r) + fldoU, (66) 

r— r-t-w 

with do = d/dr denoting an imaginary time derivative 
and Leon is a Lagrange multiplier term imposing the two 
independent U{V) gauge constraints, 

Leon = - ao{r)[erjdi(3j{T) - ttN,] 

r— r— w 

+ - /5o(r)[e„a.a,(r)-^/t^/J. (67) 

r— r-f-w 

Here we have introduced two Chern-Simons scalar poten- 
tials as Lagrange multipliers, denoted ao(i') and /?o(r), 
which live on the dual and original lattice sites respec- 
tively. 

Upon introducing another scalar potential, ao(r), living 
on the sites of the dual lattice, and collecting together 
the longitudinal and transverse parts of Oj and Cj, the 
full Euclidean action can be compactly cast into a simple 
form. In order to make the vortex physics more explicit 
we choose to re-introduce the vortex phase-field 6 within 
the Lagrangian formulation. Specifically, we shift — > 
+ cl^j^O with /i = 0,x,?/, and then integrate over both 
and 9. In this way we arrive at the final form for the 
full Euclidean Lagrangian: 

S^So + Sf + Ses + Sa. (68) 

The charge sector action Sc = '^c can be expressed 

in terms of the Lagrangian density, 

Cc^Ca + - — {do9 - ao + ao)'^ + Lkin, (69) 

with the uq — > limit enslaving oq = do9 + ao and 

U' A ■ v"^ 
= Tr(ej -y-iej{dQaj-djaa)+-^{eijdiaj-npof. 

Z TT Ally 

m 

The vortex kinetic energy terms Ckin are given explicitly 
by Tikin in Eqs. H45I47|I except with a* aj. 
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The fermion action Sf — Cf is given by, 

£f=fU^o-^|3o)fra+nf, (71) 

with Hf given in Eq. H48|) . The two sectors are coupled 
together by the electric field in the electron hopping term, 
and by the Chern-Simons action Scs — J2r=r+w 
with 



1. 



-Po{r)eijdiaj{r) 



+ -/?i(r)ey [djao{r + Xi) - doaj{r + Xi)] .(72) 



Notice that in the absence of the electron hopping term 
(te — 0) , the "electric field" ej enters quadratically in the 
action, and can then be integrated out to give Ca —^Ca, 
with 

Ca = ^(^0%- - djaof + vl{eijdiaj - npoY. (73) 

Finally, it is useful in some circumstances to treat the 
external gauge field by making the shift ^ + Ai/ir, 
which removes all coupling of Ai to the fermions, and 
furthermore leaves Ai linearly coupled to a charge "3- 
current" of the usual form, Sa = X^r with 

CA^iA^ir)Mv). (74) 

Here is the charge 3-current given explicitly by, 

Jo{r) ^ -eijOiUjir), (75) 

TT 



J,;(r) = -eij[djao(r + Xi) ~ doaj{r + Xi)], (76) 

TT 

with r = r + w. Notice that the 3-current is conserved 
as required: doJo{r) + diJiir — Xi) = 0. This form is 
useful for a variety of calculations, particularly within 
the purely bosonic RL model discussed in Sees. IIIIIIVI 
but less so in some RFL calculations best done in the 
"electron gauge" (see below) , which is incompatible with 
the above shift of e^. 

As can be seen from the equations of motion obtained 
from, 6C/Sa — and SC/6/3 ~ 0, the effect of the Chcrn- 
Simons term is to attach tt flux in a (/3) to the spinon 
(vortex) world-lines; 



(77) 



where and are the spinon and vortex three- 
currents. Here /i, v, X — 0, x, y run over the three space- 
time coordinates. 

Finally we comment on the nature of the gauge sym- 
metries of the full action S in the Lagrangian representa- 
tion. In particular, associated with the three gauge fields 
a^,afj_ and (3^ are three independent space-time gauge 
symmetries. Specifically, these are: 



2. 



Or — > Or -\- ©n 



3. 



/3m (r) 



had 



(78) 



(79) 



(80) 



with 0r, Xr and Ar arbitrary functions of space and imag- 
inary time. Because of the gauge invariance of S under 
these three distinct transformations, we are free to fix 
gauges independently for the three gauge fields. 

In addition to these three gauge symmetries, the full 
Lagrangian has two global symmetries. By construction, 
the spinon Lagrangian Cs conserves the total spin, and 
since the electrical 3-current in Eqs. I|75I76|I satisfies a 
continuity equation, 9^ — 0, the total electrical charge 
Q = Jo(r) is also conserved. 

A particularly convenient gauge choice for the gauge 
field /?p is, 



/3Kr) = -TreKr), 



(81) 



with ei defined in terms of the "electric field" ej{r) in 



Eq. (|22|l . Remarkably, in this gauge the electron creation 
operator equals the spinon creation operator. To see this, 
it is convenient to shift ao ao + ao and then integrate 
out oq, which constrains V x /? = ttV • e, or equivalently 



/3| 



Together with the above gauge choice this 



implies that Pi = —nei, so that from Eq. H55|) one has 
Cra = fra ■ We refer to this as the "electron gauge" . The 
possibility to choose a gauge within the Lagrangian for- 
mulation of the C/(l) vortex-spinon theory which makes 
fra an electron operator, suggest that it should be pos- 
sible to re-formulate the vortex-spinon Hamiltonian en- 
tirely in terms of vortices and electrons. This is indeed 
the case, as we demonstrate briefly in Appendix IdI 



III. THE ROTON LIQUID 

We first focus on the bosonic charge sector of the the- 
ory, ignoring entirely the fermions. Specifically, in the 
full Euclidean action in Eq. (|68|l we retain only the charge 
action, Sc, and the coupling to the external electromag- 
netic field, Sa- We take the fermionic density as a non- 
fluctuating constant. As discussed in Sec. Ill C II this 
constant should be determined by energetics. We as- 
sume here that the largest energy scales in the problem 
are those of the vortices, 1.6. tVf ^ Uy etc. In this case, one 
expects that vortex kinetic energy (rotonic or otherwise) 
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is minimized when the vortices experience zero average 
magnetic flux. We therefore choose the fermionic density 
equa7 to the total charge density, setting ^eijOiUj — po 
and also putting ao — 0. Note that this choice is es- 
sential to recovering an ordinary Fermi liquid state (see 
Sec. IIV AI Add. IXjl and hence is also natural in this sense. 
We remark that, while we will continue to assume the av- 
erage fermion density is equal to po in the bulk of this 
paper, we will return to another possibility - and its phys- 
ical regime of relevance - in the discussion section. 

It is also convenient to isolate the fluctuations in the 
charge density by defining, 

aj = a] + aj, (82) 

with "background" density ^eijdta'j — po- We can then 
take aj = , so that Uj — aj = dj . Of the three vortex 
kinetic energy processes which enter in Ckin, we hereafter 
and in the rest of the paper drop entirely the vortex pair 
hopping process in Eq. (|46|l putting t2v = 0. For now 
we also set the single vortex hopping processes to zero, 
putting t„ = in Eq. (|45|l . but will return to their effects 
in Section IV. Of interest here is the new roton hopping 
process, Cr = Hr in Eq. (|47|1 . which can be conveniently 
recast in the form: 

Cr = ~KrC[d] cos[A^y9 - ^{dxdy + dyUx)], (83) 

with 

C[a] = cos{e,jdihj/2). (84) 

More generally, with spinon fluctuations included one 
has C — cos[eijdi{aj — aj)/2]. Clearly, C is maximal - 
and hence the rotonic kinetic energy most negative - for 
tijdihj = 0, which is true on average for this choice of 
fermion density. 

We next choose the gauge V • a = 0, and integrate over 
ao (with uq 0). Having dropped the vortex hopping 
processes, the remaining charge Lagrangian is then given 
by, 

Cc = Cpi+Cg, (85) 

with 

^Pi = ^[(c'oaj)^ + Voi^^j^^djf], (86) 

Cg = -^{djdoOY - KrC[a\ cos[/\xvO - \{did2 + 820.1)]. 

(87) 

To analyze the phases of this model it is instructive to 
represent the Lagrangian Cg in Hamiltonian form by re- 
introducing the vortex number operator, N^: 

He = ^^^riV,,y(r-r') 

r,r' 

-Kr^C[a]cos(Aj;j^0r - '^{dxO'i +dydx)). (88) 



A 

T 




FIG. 1: Schematic phase diagram of the Hamiltonian in 
the charge sector, He = Hg + Hpi, with zero vortex hop- 
ping strength, — 0. When the roton hopping vanishes, 
Kr = 0, He describes a classical logarithmically interacting 
vortex gas, and has a superconductor-to-normal transition 
at a Kosterhtz-Thouless temperature, Tkt ~ Uv- With in- 
creasing roton hopping, Tkt decreases being driven to zero 
at K* ~ Uv where there is a quantum phase transition from 
the superconductor into the Roton liquid phase. 



Again V{r) is Fourier transform of the discrete inverse 
Laplacian operator, V(k) = 1//C^(k), with /C^(k) = 
2(1- cos /c,). 

The first term in Eq. (|88|l describes a logarithmically 
interacting gas of (integer strength) vortices moving on 
the dual 2d square lattice. When = 0, this model 
will undergo a finite temperature Kosterlitz-Thouless 
transitions^ from a high temperature "vortex plasma" 
into the low temperature "vortex dielectric" . This cor- 
responds, of course, to a transition into a superconduct- 
ing phase. With = the Kosterlitz-Thouless tran- 
sition temperature will be set by the vortex interaction 
strength, Uv But upon increasing the strength of the ro- 
ton hopping, one expects the transition temperature to 
be suppressed, and for ^ Uy to be driven all the way 
to zero. Thus, at zero temperature, upon increasing the 
single dimensionless ratio, K^/uy, one expects a quan- 
tum phase transition out of the superconducting phase 
and into a new phase (see Fig. 1) - the "Roton Liquid". 



A. Harmonic theory and Excitations 

To access the properties of the Roton Liquid (RL), we 
consider k,. ^ Uy , where it is presumably valid to expand 
the cosine terms in Eq. (|88|l for small argument, giving 
Hg = Hrot + with 

r,r' r 
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K . '1 



(89) 



With this expansion, it is no longer legitimate to restrict 
6 to the interval [0, 27r]. Consistency then dictates that 
the eigcnvahics of the vortex number operator no longer 
be restricted to integers, but allowed to take on any real 
value from [—00,00]. 

The full Roton Liquid Hamiltonian. Hrl = Hpi +-ffrot 
is quadratic and can be readily diagonalized. This is 
most conveniently done by returning to the Lagrangian 
framework, described now by 



1 

2uv 



i^xyO- Tiid^ay + dya^)]'', (90) 



with vo ~ ^1^0+ KrUy/4:. To proceed to describe the 
normal modes of this quadratic Lagrangian, we define 
Fourier transforms 



0(r,r) = ^ 
0(r,r) = ^ 



_ik-r— ioJnT 



0(k,w„), 



0(k,a;„ 



(91) 
(92) 



for fields O, on the original and dual lattices, respec- 
tively. Here integration = J d^k/(27r)^ is taken over 
the Brillouin zone Ifel < TT and = /^cL;„/(27r) 
defines the integration measure (at zero temperature) for 
the Matsubara frequency ojn- At non-zero temperature, 
one simply replaces l3~^ Sa;„' ^^^^ = 27rn//3. 

It is moreover convenient to define 



/C,(k) = -i(e^'=^-l), 



(93) 



so that upon Fourier transformation, the discrete deriva- 
tives behave intuitively. 



i/Cj(k), 



(94) 



and of course do —ioOn as usual. We also introduce the 
transverse gauge field: 



ai(k) 



/C(k) 



(95) 



with /C2(k) = J2j l^j(k)|^ and |/Cj(k)| = 2| sin(A:j/2)|. 

To now diagonalize jCrl, it is convenient to define a 
real two-component field via 



a(k,a;„) = 

9{\s.,U}n) = 



fu 



'Ti(k, Wn), 



/C(k) 



^T2(k,a;„). 



(96) 
(97) 



Then the action, Srl = /r 

Srl = \ /T„(k,w„)C;-i(k,a;„)T^(-k,-a;„),(98) 



with 



G^6. 



a/3 



a/3 



a/3 



a/3 



(99) 



where a is the usual vector of Pauli matrices. Here we 
have defined, 



G' 



2 t^2 , 2 \l^xlCy\'^ 



f(i/c.r-i/c,n^ 



(100) 
(101) 
(102) 



with vi = sjKrUv and /Cj = 2sin(A;j/2). The poles in 
Gafi at = itjjn = ±a;p;,±a;rot describe two types of 
collective modes. 

The first excitation is a plasmon with a renormalized 
dispersion, 



with velocities, 



(103) 



v± = 



±v2/4. 



(104) 



The plasmon frequency vanishes at the center of the Bril- 
louin zone, k = 0, and in the absence of long-ranged 



Coulomb interactions disperses linearly u) 



pi 



Vpl 



Ikl at 



small wavevectors, kj 0. But the associated plasmon 
velocity, Vpi{(j)), depends upon the ratio, ky/kx = tan((/)). 
In particular, along the zone diagonals with ky = ±kx 
the velocity is minimal and unaffected by the vortices 
with Vpi = ^0, whereas it takes it's maximum value, v+, 
along the kx or ky axes. 

This upward shift in the plasmon frequency is due to 
a "level repulsion" with the second collective mode - the 
gapless roton, which disperses as. 



'^rot(k) = 



1 



- \/vi}c^ + vlvli\K:x\^-\}Cy\^y 



(105) 

For \kx\ <^ 1 and fixed ky the roton dispersion vanishes, 

l^rot ^ Vrot\kx\, with 



Vrot = 1^1- 



(106) 



Remarkably, the Roton Liquid phase supports a gapless 

"Bose surface" of roton excitations, along the kx = and 
ky = axes. These roton excitations describe gapless and 
transverse current fluctuations, which are obviously not 
present in a conventional bosonic superfluid. 

With long-range Coulomb interactions present one 
would have simply, Vq — > i)o(k) giving the familiar 

2d plasmon dispersion, ujpi ^ addition, the ro- 

ton velocity Vrot becomes dependent upon ky. We note 
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in passing that the roton velocity is in either case deter- 
mined not only from the dynamics of 9 but also that of 
a, as is evident from its dependence upon i;o/w+. It will 
be sometimes instructive in the following to consider the 
simple limit ~ "Do OO; in which the spatial fluctua- 
tions of a vanish and the roton mode is entirely decoupled 
from a. 



pair field, Br,r+x, , 

G?J(r,r) = (So,i.(0)sL+x,W)' (109) 

behaves similarly, and will be discussed at the end of this 
sub-section. 



B. No Meissner effect in roton liquid 

We now employ the Gaussian theory to examine some 
of the electrical properties of the Roton liquid phase. 
Consider first the response of the RL phase to an ap- 
plied magnetic field. In the presence of a magnetic field, 
B = eijdiAj, there is an additional term that one must 
add to the Lagrangian, which from Eq. ifT^ takes the 
form: 



Ca — —aoB. 

IT 



(107) 



If oo is integrated out from Cc in Eq. (|69|l with this 
additional term present, the Hamiltonian, Hg{Nr,9r) in 
Eq. becomes simply, He{Nr ~ Q,). As expected, 
the vortex density will become non-zero in the presence 
of the magnetic field. Since the vortex number operator 
in Hq has integer eigenvalues, it is not generally possible 
to shift away the applied B-field. But in the Roton Liquid 
phase (at zero temperature) where the cosine term can 
be expanded to quadratic order (as in Hrot): the vortex 
number operator has a continuous spectra, and one can 
formally eliminate the B-field by shifting, N^ + ^B, 

for all r. Since the ground state energy of the RL phase 
is thus independent of the applied B-field, both the mag- 
netization and the magnetic susceptibility, x = dM/dB 
vanish. Unlike in a superconductor, where X = "j^i 
there is no Meissner effect in the Roton Liquid (strictly 
speaking, there is never a Meissner effect in a single two 
dimensional layer, but one can consider an infinite stack 
of electrically decoupled but magnetically coupled lay- 
ers, which would exhibit a Meissner effect when the lay- 
ers are true 2d superconductors, but not when they are 
RLs). Physically, since the RL phase supports gapless 
roton excitations, the state cannot screen out an applied 
magnetic field. 



C. Off-Diagonal Quasi-Long- Range Order 

We next consider the Cooper pair propagator in the 
Roton Liquid, 



G'^ir, ~ r2,Ti - T2) = (i?r,(TlX(r2)), 



(108) 



with the Cooper pair destruction operator B^ given in 
Eq. H37|l as an infinite product of exponentials, e^'^^'^J , 
running along the string. The propagator for the d-wave 



1. Equal time correlator 

We consider at first the equal time correlator, with 
Ti = T2 . The path independence of the string rests on the 
condition, (V-e)(r) = Nr, with mtet^er vortex number, N^. 
Unfortunately, within the tractable harmonic approxima- 
tion valid for most quantities in the roton liquid phase 
(with cosine terms in the roton hopping expanded to 
quadratic order), the condition of integer vortex number 
is not satisfied, and the results for G"^^'(r, 0) depend upon 
the choice of string. We believe that the correct behavior 
can be extracted by taking the string running along the 
straightest and "shortest" (using the "city block" met- 
ric \xi — X2\ + \yi — 2/2!) path between the two points, 
ri and r2. As we shall see, the Cooper pair propagator 
calculated in this way has an anisotropic spatial power- 
law decay. Preliminary calculations suggest that, once 
perturbative corrections to the harmonic approximation 
(using the formalism established in Sec. II V B^ are taken 
into account (even if they are irrelevant in the renormal- 
ization group sense), simple variations in the string do 
not modify the power-law decay of G"^P(r,0), but only 
change the (non-universal) prefactoriii 

We take ri — r2 = Xx + Yy, and r2 = w, with in- 
teger X,Y > 0. Then, upon expressing the correlator 
G'^P{X, Y, 0) as an imaginary time path integral, one ob- 
tains an extra term in the Euclidean action, S = J ^ 
in Eq. with 

C + iej{r,T)Jj{r,T). (110) 
The c-number "source" field is given as, 

X-l Y-l 
Jji^^'^) = 2TrS{T)[dj2 ^ Sx,x' + lSy,0- Sjl ^ Sx,xSy,y']. 

x'=0 y'=0 

(111) 

The Fourier transform is simply, 



J,-(k,Lj„) = 27r[(5j2 



ilCx 



iJCy 



(112) 

Integrating out the electric field ej in the gauge W ■ a — 
and decomposing the source fields into transverse and 
longitudinal parts, J7t — ieijK-iJ'j / IC, Ji = iK*Ji/K, one 
obtains the result 

G^P(A,r,0) = (exp{- / [^(j-^Ti- J,T2) 
1 
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where the Gaussian average over T is to be taken with 
respect to Srl in Eq. Performing this Gaussian 

integral, one obtains 



where 



G^p(x,y,o) = cxp[-ri-rt-r,t], 



\Jl\\l-ulG22) ^ 

\Jt\Hl-LolG,^) 



Tit 



(114) 

(115) 
(116) 
(117) 



with Gij given in Eq. 

Investigation of Tt and Yu shows that the correspond- 
ing integrands are non-singular at small or /c^, and 
hence go to finite limits for large |X| and/or large 
They will thus affect only the amplitude of the Cooper 
pair propagator at large distances, and we henceforth ne- 
glect them. Singular behavior at long distances does arise 
from F; , in line with the intuition that it is vortex fluctu- 
ations which disrupt the superconducting phase, since Ji 
couples to the longitudinal electric field, which through 
V • e = describes the vorticity. To evaluate F; , we first 
perform the frequency integration to obtain 



, i)oi.i|/C./C,| + ff|/C./C,|V/C^ 



(118) 



Next, we explicitly express the square of the longitudinal 
string. 



/C2 

-t-2Re{(l 



1 



JC. 



)(1 



ll-e 



(119) 



The first two terms in Eq. (|119(l are singular for small ky, 
kx respectively for very large Y,X, leading to a logarith- 
mic dependence when inserted in Eq. (|118|l . The final 
term in Eq. (|119|l . by contrast, is singular only for both 
kx,ky small, and this singularity, inserted into Eq. (|118|1 . 
is weak and integrable. Extracting the logarithmic parts, 
one finds 



F,^A,(ln|X|+ln|y|), 
for \Xl \Y\ > 1, with 



2. Unequal time correlator 

We now consider the Cooper pair propagator at un- 
equal times. Unfortunately, it is difhcult to produce 
a simple and general calculation for arbitrary spatial 
and time separations. In particular, clearly, by square 
symmetry, we expect G^P(X,y,T) = G''p{Y,X,t). Any 
choice of strings, however, necessarily creates an asym- 
metry between the two spatial directions. As we have 
been unable to resolve this dilemma, we instead focus on 
the simple case in which the pair is created and annihi- 
lated on a single row of the lattice, i.e. G'^'''{X, 0, r). We 
will see that this correlator decays as a power law both 
in space and time. 

To proceed, we take ri — r2 = Xx, r2 — w, ti — r, 
T2 = 0, with X,Y > 0. With this choice, the string in 
Fourier space becomes 

2 _ g-ikxX + lLJnT 

Jj{k,UJn) = 27TSj2 . (123) 

Repeating the same manipulations as above, one again 
obtains (with negligible contributions from the transverse 
part of the string) 



G^f(X,0,r)^exphf,], 



with 



A, 
2 



r,^^Hx^ + vl,r% 



(124) 



(125) 



where for simplicity we have taken Vrot independent of ky 
(for a non-trivial v{ky), the logarithm is simply averaged 
uniformly over the ky axis). This gives 



GciX,0,T) 



const. 



(X2+«2^,r2)Ac/2- 



(126) 



Note that this power-law form implies a power-law local 
tunneling density of states for Cooper pairs, Pcpi^) ^ 



We conjecture that the full correlator satisfies a simple 
scaling form with "z=l" scaling: 



G^P^X, Y, t) 



const X 



Y 



IXlAclylA, 



VrotT VrotT 



). (127) 



(120) Combining our two calculations above implies 5(x, 0) 

(xV(x' + i))^=/'- 



(121) 



3. ODQLRO of the d-wave pair field 



Hence we have 



G''P{X,Y,0) 



const 



Xl^^lYl 



(122) 



This establishes that the Roton liquid phase has off- 
diagonal quasi-long-ranged order (ODQLRO) at zero 
temperature. 



We now briefly discuss the analogous ODQLRO of the 
d-wave pair field described by G^j'(r,T) in Eq. (|l()9|l . 
This quantity is plagued by the same string ambigui- 
ties as the local Cooper pair propagator, but to a larger 
degree, since it involves two separate strings emanating 
from the two sites shared by the initial bond and end- 
ing at the two sites shared by the final bond. Although 
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we are confident G^j'(r, r) has a power-law form consis- 
tent with ODQLRO, we are unable to determine the pre- 
cise nature of these correlations with reliability. For in- 
stance, consider the equal time pair field correlator for 
two bonds along the x-axis, Gj^(r,0). For two bonds in 
the same row, r = Xxi, the strings can be chosen to 
line all in the same row, and since the logarithmic di- 
vergence controlling the ODQLRO arises from small 
in this case, the power law exponent is unchanged, i.e. 
Gif(X,0,0) ~ Const/|X|^-, with given above. We 
believe that, since this choice of string is by far the most 
natural, this is probably the correct result. If, instead, 
we choose to separate the two pair fields along a single 
column, r = Fx2, then the two strings involved cannot 
be taken entirely atop one another. Different choices for 
the strings then give different results. For instance, mak- 
ing the symmetric choice of two parallel strings (each of 
strength tt rather than 27r) gives a decay exponent re- 
duced from Ac to Ac/2 in the Y direction, while choos- 
ing the strings to overlap everywhere except the two ends 
reproduces the previous exponent A^ without any reduc- 
tion. Since we are unable to reliably resolve this am- 
biguity, we are unable to determine the exact form of 
the d-wave pair field correlator. Instead, we will take the 
pragmatic approach of approximating the correlations by 
those of the local pair field, G^j'(r,T) « G'=P(r,r). It 
should be understood that the decay exponent Ac may 
need to be renormalized and/or the correlator corrected 
slightly to obtain detailed results for a specific model. 



with 

n?,(k,-„)-5^, (130) 

as |k| — > 0. The dependence of 11^^- on the orientation of 
k is due to the fact that we have assumed i? = 0, which 
according to Faraday's law requires VxE = — i9tB = 
0. Hence we must choose k parallel to the electric field 
A — E/(— icij). Thus we extract an isotropic conductivity 
tensor cr" — Sija'^, with 



characteristic of a system with no dissipation. 

The above conclusion for the quadratic RL Lagrangian 
is, however, modified by the vortex hopping terms. As 
we detail in the next section, the effects of a small vor- 
tex hopping term depend sensitively on the parameters 
that enter in the harmonic theory of the roton liquid - 
in particular the dimensionless ratio Uy/ k.^- There are 
two regimes. When this ratio is larger than a critical 
value, vortex hopping is "relevant" and grows at low 
energies destabilizing the roton liquid phase. On the 
other hand, for small enough u^/nr the vortex hopping 
strength scales to zero and the roton liquid phase is sta- 
ble. In this latter case, the effects of vortex hopping on 
physical quantities can be treated perturbatively. In par- 
ticular, we find that the conductivity in the roton liquid 
diverges as a power law in the low frequency and low 
temperature limit. 



D. Conductivity in the harmonic theory 

Given the above result of ODQLRO, it is natural to 
expect a very large and perhaps infinite conductivity in 
the Roton liquid. Indeed, neglecting the effects of vor- 
tex hopping, Cy in Eq. (|^ . one can readily see (e.g. 
from Eq. (|88|l l that the total vortex number on each row 
and each column of the 2d lattice is separately conserved. 
Thus it is impossible to set up a vortex fiow, and hence 
by the Josephson relation to generate an electric field. 
Using the quadratic roton liquid Lagrangian Lrl, this 
expectation can be directly confirmed. In particular, we 
can integrate out all dynamical fields (0, a) from the La- 
grangian leaving only the external gauge field , and 
thereby extract the polarization tensor and conductiv- 
ity. This is most conveniently carried out in the gauge 
Aq = 0, and assuming no magnetic field dxAy—dyAx = 0. 
In this case, one may write 

- I '^^"^^ C^) e-''--Wa(k, u^)A, (-k, -^.„). 
■/k,i^„ 7r/C(k) 

(128) 

Integrating out a and 9 using the Green's function Gn in 
Eq. (|99|l , one obtains, in the limit |k| 0, the effective 
action 

Sa = \ I n^,(k,c^„)A,(k,c^„)A,(-k,-c^„), (129) 



IV. INSTABILITIES OF THE ROTON LIQUID 

We first consider the instabilities of the roton liquid 
due to the presence of a vortex hopping term, and exam- 
ine the effects of such processes on the electrical trans- 
port. In the subsequent subsection we consider the le- 
gitimacy of the harmonic expansion required to obtain 
the quadratic RL Lagrangian. This is achieved by per- 
forming a "plaquette duality" transformation^ «, where it 
is possible to address this issue perturbatively. Again 
we find two regimes depending on the parameters in the 
quadratic roton liquid Lagrangian. A stable regime for 
small Uy/Kr wherein the harmonic roton liquid descrip- 
tion is valid, and an instability towards a superconduct- 
ing phase when this ratio is large. 

A. Vortex hopping 

The analysis of the stability of the RL to vortex hop- 
ping is mathematically nearly identical to the stabi lity 
analysis of the Exciton Bose Liquid (EBL) of Ref. [IJ 
with respect to boson hopping. Taking over those meth- 
ods, we note that the vortex hopping operator exhibits 
one-dimensional power-law correlations, 

{e'9y9m^-^dye,^Ar))^^ = <5,,o7^(x, r), (132) 
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with a power law form at large space-time separations, 



7^(a;,T) 



(133) 



where r = xx + yy and c is a dimensionless constant. 
Here, TZ{x,t) is essentially the single roton Greens func- 
tion, describing the space-time propagation of a roton 
with a dipole oriented along the y axis. When calculated 
using the RL Lagrangian, one finds 

(134) 

Simple calculations show that this power-law behavior is 
not modified by including the fiuctuating ay field in the 
vortex hopping operator, Tj(r,T) — e^^^^^'^"'''^''>\ i.e. 



{fyir,0)f;{r + r,T))nL 



Sy,0 



(a;2 



(135) 



with only a change in the prefactor. Notice that the 
exponent A„ characterizing the power law decay of the 
roton propagator is inversely proportional to the anal- 
ogous exponent Ac in Eq. (|121|l which gives the power 
law decay of the Cooper pair propagator. Indeed, for the 
Roton liquid Lagrangian studied here we find the simple 
identity, Ay Ac = 1/2. But with inclusion of other terms 
in the original Hamiltonian such as the spinous or fur- 
ther neighbor roton hopping terms, this equality will be 
modified. 

The arguments of Ref. imply that the vortex hop- 
ping term is then relevant for < 2. In this regime, 
the vortex hopping strength grows large when scaling to 
low energies, and one expects the vortices to condense 
at zero temperature. In this case it is legitimate to ex- 
pand the cosine term in Eq. H45|) and one generates a 
"dual Meissner efi^ect" where the gauge fields that are 
minimally coupled to the vortices become massive. In 
the presence of spinous this leads to a mass term of the 
form, Cy ~ ^(oj — QijO^j which confines one unit of elec- 
trical charge to each spinon presumably driving one into 
a Fermi liquid phase. If we ignore fluctuations in the 
spinon density, or drop the spinous entirely retaining a 
theory of Cooper pairs, the resulting phase is a charge 
ordered bosonic insulator. The nature of the charge or- 
dering will in this case depend sensitively on the com- 
mensurability of the Cooper pair density (po/2) with the 
underlying square lattice. In the simplest commensurate 
case with one Cooper pair per site {po = 2), a featureless 
Mott insulating state obtains. 



1. Electrical Resistance in the roton liquid 

When A^ > 2, on the other hand, vortex hopping is 
"irrelevant", and the effects of the hopping on physical 
quantities can be treated perturbatively. Despite its irrel- 
evance, we expect the vortex hopping to strongly modify 



the Gaussian result for the conductivity by introducing 
dissipation. To understand how this occurs, it is instruc- 
tive to first consider the simple limit alluded to earlier 
in which vq oo. In this limit, the longitudinal den- 
sity fluctuations described by d at non-zero wavevector 
are suppressed, and the roton mode is purely captured 
by the 9 fleld. The zero wavevector (but non-zero fre- 
quency) piece of a, however, remains non-zero in this 
limit and can be used to calculate the conductivity in an 
RPA fashion. In particular, we take into account the ro- 
ton fluctuations and their associated dissipation induced 
by vortex hopping by calculating the effective action for 
a, A upon integrating out 6 to second order (the lowest 
non-trivial contribution) in ty. Starting then with the 
Lagrangian, Crl + Ca + ^v, expanding to second order 
in the vortex hopping action Sy , and integrating over the 
9 field and the gauge field a(k ^ 0) with vq ^ oo gives, 

S^A = E / [^(^o'^^)' - ^6,,a,(r)9oA,(r - xj)] 

r— r+w ^ 



, o(2) 



(136) 



where 



^2 r 

Sa^A = - Y E / , (C0s(az6l - ai),r C0s{dj9 - aj )r'r')e ' 
r.r' "^'"^ 

(137) 

where {■)g indicates the average with respect to the Gaus- 
sian action for 9. From Eq. H132|) . one can carry out this 
average to obtain. 



-(2) 
'a, A 



X cos[ay(r, r) — dy{r + xx, r')] + {x ^ y) 



438) 



with the roton propagator TZ{x,t) given as in Eq. (|133|l . 
except with A„ — > Uy / / {^tt) in this limit. Follow- 
ing the usual RPA strategy, we expand S*'^' to quadratic 
order in a to obtain, 



:(2) 



4 ./k.. 



'R{uJn)\dj{uJn)\ , 



(139) 



with the definition, TZiujn) — 7?.(0) — 7?.(cj„) and, 

7^(cJ„) =7^(A:, = 0,w„)- (140) 

The k = limit is valid when — > cxd . At low frequencies 
one has, 



' Vrot sin -1(1 -f 7) 



(141) 



with > a dimensionless constant. Here we have 
retained explicitly the leading singular frequency depen- 
dence and dropped analytic terms consisting of even pow- 
ers of LOn- The exponent 7 is defined as. 



7 = 2A^ 



(142) 
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and in the stable regime of the RL phase, 7 > 1. 

Finally, upon integrating out dj, one obtains a renor- 
malized electromagnetic response tensor 



(143) 



It is now straightforward to extract the conductivity 
by analytic continuation. 



a(w) 



One obtains an appealing Drudc form: 

1 



-iuj{TT'^/Uy) + i^^7?.,.et(w)/w 



(144) 



(145) 



with the retarded propagator obtained from analytic con- 
tinuation: 



7^ret(t^) = 7?.(t^n)| 



(146) 



The non-analytic frequency dependence oi TZret{^) con- 
tributes to the dissipative (real) part of the resistance 
(per square), R{uj) = Recr~^(ct;), which is quadratic in 
the vortex hopping amplitude, if,: 



R{lo) 



2 7-2 



2u) 



-ImTZreti'-^) = 



Vrot 



(147) 



Note that at the point for which vortex hopping is just 
marginal, 7=1, the resistance becomes linear in fre- 
quency Re(T(w) ^ 1/uj. 

We can readily extend this result to finite tempera- 
tures, by using the finite temperature roton propagator: 



7?.(t) = c-y 



1^1 Vrot 13 

sin{-KT / j3) 



1+2 



(148) 



The retarded roton propagator follows upon analytic con- 
tinuation, and can be most readily extracted by using the 
identity: 



/TO7^re^(w) = sinh(/3w/2) / dte*"*7^(r 



it). 
(149) 

Upon combining Eq. (|147I149|I we thereby obtain a gen- 
eral expression for the finite temperature and frequency 
(dissipative) resistance in the roton liquid phase: 



R{lu,T) 



(vrt„)2 



ttT 

Vrot 



with a universal crossover scaling function, 
27 



R^{X) 



r(2 + 7) 



^ 2 X 



(150) 



(151) 



interpolating between the d.c. resistance at finite tem- 
perature and the T = a.c. behavior. Since R^{X 0) 



is finite, the d.c. resistance varies as a power law in tem- 
perature: R{T) ^ T'' . At the boundary of the RL phase 
with 7 = 1 , a linear temperature dependence is predicted. 
At large argument. 



i?-y (yX 



r(2 + 7) 



(152) 



so that the resistance crosses over smoothly to the zero 
temperature form, R{uj,T = 0) ~ |a;|'>'. 

This RPA treatment has the appeal that it produces 
the natural physical result that the effect of the weak (ir- 
relevant) vortex hopping is to generate a small resistivity 
^ t1 . Formally, it is correct for vq = 00 because the RPA 
reproduces the exact perturbative result for the electro- 
magnetic response tensor to O(t^) in this case. Unfortu- 
nately, when the spatial fluctuations of hj are not negligi- 
ble, i.e. for Do < 00, even the 0{t^) term is not obtained 
correctly. More generally, the fluctuations of hj and 9 
must be treated on the same footing. Therefore in the 
general case we instead integrate out hoth fields and ob- 
tain more directly the correction to litj to 0{if). The 
calculations are described in appendix IeI. This does not 
yield the appealing "Drude" form in Eq. 1145|) but in- 
stead the Taylor expansion of Eq. (|143|l to O(t^), 



n 



(2) 



2, ,2 



tin 



2v: 



2A„-1 
rot 



|2A„-4 kjkj 



(153) 



except with the scaling dimension A„ , given explicitly in 
Eq. H134|) . now fully renormalized by the plasmon fluc- 
tuations. Provided the vortex hopping is irrelevant, this 
is sufficient to recover properly the low-frequency behav- 
ior of the resistivity, Eq. (|147|l to 0{tl). In particular, 
formally inverting the perturbative result for a{Lu,T) to 
O(i^), we infer the appropriate d.c. dissipative resistance 
i?(T) ~ tlT-y, with 7 = 2A„ - 3. 



B. "Charge Hopping" 

In this subsection we examine the legitimacy of the 
"spin wave" expansion employed in Sec. Ill A to obtain 
the Gaussian Roton liquid Lagrangian. This is most 
readily achieved by passing to a dual representation, 
which exchanges vortex operators for new "charge" op- 
erators. This procedure is a quantum analog of the map- 
ping from the classical 2d ATF-model to a sine-Gordon 
representationi&, the latter suited to examine correc- 
tions to the low temperature "spin wave" expansion. As 
we shall find, there are parameter regimes where the 
"charge" hopping terms are irrelevant and the RL phase 
is stable. But outside of these regimes, the "charge" 
quasiparticles become mobile at low energies and con- 
dense - driving an instability into a conventional super- 
conducting phase. Throughout this subsection we will 
drop the vortex hopping term, Hy, focusing on the pa- 
rameter regimes of the Roton Liquid phase where it is 
irrelevant (i.e. A„ > 2). 
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1. Plaquette Duality 

To this end we now employ the "plaquette duality" 
transformation, originally introduced in Ref. [T2| in the 
context of the Exciton-Bose-Liquid phase. Consider the 
charge sector of the theory, with Hamiltonian He = Hpi + 
He + H^. This Hamiltonian is a function of the vortex 
phase field and number operators, 9^, (living on the 
sites of the dual lattice), as well as the (transverse) gauge 
field a* and it's conjugate transverse "electric field", e*. 
The plaquette duality exchanges Or and Nr for a new 
set of canonically conjugate fields which live on the sites 
of the original 2d square lattice. The two new fields, 
denoted 4>r and hr, are defined via the relations, 

7riVr+w = ^xy(i>r, (154) 
Trflr = Aa;y0r-w (155) 

Although ipr and fij. are conjugate fields satisfying, 

[rir,(^r'] = «<^rr', (156) 

they can not strictly be interpreted as phase and number 
operators since the eigenvalues of 4>r = ttto for arbitrary 
integer m, whereas Trfir is 27r periodic. It is important 
that (pr and rir not be confused with the "chargon" phase 
and number operators introduced in Section II which 
were denoted (pr^Tir - without the tildes. As we discuss 
below, e*"^' does in fact create a charge-like excitation, 
but it is not the chargon. 

Under the change of variables, Ho{9, N) —f H^{<f), h), 

H^ = Hu- Kr ^ cos[7rnr+w - ^{dxay + dyCi^)], (157) 

r 

with the definition, 

Hu^^Y. ^xy^rA,y4,,,V{r - r'). (158) 

In the Roton Liquid phase the cosine term in i?^ is 
expanded to quadratic order, and H^ + Hpi Hrl. To 
be consistent, both fi and (j) must then be allowed to 
take on any real value, and it is convenient to pass to a 
Lagrangian written just in terms of (p: 

T - T ^H , 1 vlMil! 

Lrl - Lpi + Hu + -2_^ 

r ^ 

+ ^y^do4>r+^{dxay + dyar,). (159) 

r 

In this dual form the Roton Liquid Lagrangian depends 
(quadratically) on the field 0r, which lives on the sites of 
the direct lattice, and the (transverse) gauge field, aj(r), 
defined on the links of the dual lattice. 



2. Superconducting Instabilities 

This dual formulation is ideal for studying the legiti- 
macy of the "spin- wave" expansion needed to obtain the 
quadratic RL Lagrangian. The crucial effect of the spin- 
wave expansion was in softening the integer constraint on 
the eigenvalues of ^/tt, allowing (p to take on all real val- 
ues in LuL. It is, however, possible to mimic the effects 
of this constraint by adding a potential term to L^l of 
the form, 

La = -A^cos(20r)- (160) 

r 

When A — !■ cx) the integer constraint is enforced, whereas 
the RL phase corresponds to A = 0. Stability of the 
RL phase can be studied by treating A as a small per- 
turbation to the quadratic RL Lagrangian. But as dis- 
cussed in Rcf. ;T3, one should also consider other per- 
turbations to Lfjj^ which might be even more relevant. 
Generally, one can add any local operator involving 0r 
at a set of nearby spatial points which is 27r periodic in 
20, and satisfies all the discrete lattice symmetries (i.e.. 
translations, rotations and parity). For example, terms 
of the form cos{2£(f>) for arbitrary integer £ are allowed, 
although these will generically become less relevant with 
increasing £. as we shall see, for our "minimal" model of 
the RL phase, the most relevant perturbation is of the 
form, 

Lt = -tcJ2Yl cos(29,0r). (161) 

Before studying the perturbative stability to such oper- 
ators, we try to get some physical intuition for the mean- 
ing of the operator e**^. From the commutation relations 
in Eq. H156(l . the operator e"^"' increase fir by one, and 
creates some sort of quasiparticle excitation on the spa- 
tial site r. Since the perturbation in Eq. H160() changes 
the number hr by ±2, the total number of these quasipar- 
ticles, htot = "-^^ conserved, but the "complex 

charge", Qc ~ e*'^"*°* is conserved. The perturbation in 
Eq. H161() can then be interpreted as a "charge hopping" 
process. To get some feel for the nature of the quasipar- 
ticle, it is instructive to introduce an external magnetic 
field, B — eijdiAj, which enters into H^ above via the 
substitution, A^yipr Axy4>r — Br+w A spatially uni- 
form field, B, can readily be removed from the Gaussian 
RL Lagrangian by letting 

0r -> 0r + Bxy. (162) 

But then B appears in the cosine perturbations, in "al- 
most" a minimal coupling form. For example, one has 
the combination dx4> — 2Ax where we have chosen the 
gauge Ax — —By/2 and Ay = Bx/2, suggesting that 

e"'^ carries the Cooper pair charge. But the y— derivative 
enters as dy(f> + 2Ay - with the wrong sign. Thus, this 
quasiparticle is not a conventional electrically charged 



19 



particle. Nevertheless, as we show below, condensation 
of the quasiparticle with (e*'^) 7^ does drive the RL 
phase into a superconducting state. 

To evaluate the relevance of the various perturbations 
in the RL phase, requires diagonalizing the associated 
action, Srl- In momentum space one has, 

(163) 



lly I K^x I 



Daa [<^„ 

U„ 



D. 



acj) 



c.„(/Cg-/Cg)e"-- 
27r/C 



(164) 



(165) 



(166) 



with kfj_ — (k, cjn). Evaluation of the two-point function 
of e^"^ then gives. 



with scaling dimension, 




(167) 



(168) 



We note in passing that this power-law form of the 
"charge" operator (e^**^) two-point function in Eq. p67|) 
differs from exp(— In^ r) behavior of the corresponding 
object in the Exciton Bose liquid of Ref. IT^ due to 
the long-range logarithmic interactions between vortices. 
Stability of the RL phase requires Ac > 1. Evaluation of 
the two-point function for the "charge" hopping opera- 
tor, e^'^""^, gives. 



^gi2a^ 0, (r) g-i2a„ 00 (0) ^ 



(a;2 



(169) 



with r — a;x + j/y, and with the same scaling dimen- 
sion Ac as above. However, since this two-point function 
decays algebraically in two (rather than one) space-time 
dimensions, the perturbation tc is relevant for Ac < 2. 

When Ac < 2, the "charge" hopping process will grow 
at low energies, and will destabilize the roton liquid 
phase. Not surprisingly, the resulting quantum ground 
state is superconducting. Indeed, the exponent Ac in 
Eq. lfTB5|) above is in fact identical to the exponent char- 
acterizing the power law decay of the Cooper pair prop- 
agator in Eq. (|121|l . so that for small Ac the Roton liq- 
uid phase is already "almost" superconducting. More- 
over, when the vortex core energy greatly exceeds the 
roton hopping strength, ^ Kr, the Hamiltonian Hg in 
Eq. (|S5|) is deep within it's superconducting phase. This 
limit precisely corresponds to Ac <IC 1, the limit where 



the "charge" hopping is strongly relevant. More directly, 
when the "charge" hopping strength grows large, the 
field (j) gets trapped at the minimum of the cosine poten- 
tials in Eqs. H16UI161|I . and it is legitimate to expand the 
cosine potentials to quadratic order. Once massive, the 
expectation value (e**) ^ - and the charge quasiparti- 
cle has condensed. Moreover, setting = in Eq. H161|l 
in the presence of an applied magnetic field will generate 
a term of the form. 



(170) 



indicative of a Meissner response. 

In the resulting superconducting phase, the rotons - 
gapless in the RL phase - become gapped. This follows 
upon expanding the cosine term in Eq. (|160|l . Cx = 2A^2, 
which leaves the roton liquid Lagrangian quadratic, and 
allow one to readily extract the modified roton disper- 
sion. For fcj. — !■ at fixed ky, this gives: 



w™t(k) = y W^otfc2 + ' 



(171) 



with the "roton mass gap" given by, rrirot — 27r^ k^A. 

Since the product AcA^ = i, it is not possible to have 
both the vortex hopping and the "charge" hopping terms 
simultaneously irrelevant. 



V. THE ROTON FERMI LIQUID PHASE 

We now put the fermions back into the description of 
the Roton liquid. We first consider setting the explicit 
pairing term in the fermion Hamiltonian 7i/ in Eq. I|48|) 
to zero: A = 0. As in Sec. IIIII we will assume for the 
most part (with the exception of Sec. lVIIIB3l in the dis- 
cussion) that the equilibrium fermion density is equal to 
the charge density po- This choice naturally minimizes 
Coulomb energy and vortex kinetic energy, as discussed 
therein. As we shall see, in this way we will arrive at a de- 
scription of a novel non-Fermi liquid phase - the Roton 
Fermi liquid - which supports a gapless Fermi surface 
of quasiparticles coexisting with a gapless set of roton 
modes. We then reintroduce a non-zero pairing term, 
and study the perturbative effects of A. We argue that, 
when the scaling exponent that describes the decay of the 
off-diagonal order in the Roton liquid is large enough, 
Ac > A* > 3/2, the explicit pairing term is perturba- 
tively irrelevant, and the RFL phase with a full gapless 
Fermi surface is stable. 

Even in the absence of the explicit pairing term - which 
couples the fermionic and vortex degrees of freedom in a 
highly non-linear manner - the rotons and quasiparti- 
cles interact through (three-)current-current interactions 
"mediated" by the and fields. Although these in- 
teractions are long-ranged for individual vortices, they 
are not for rotons, which carry no net vorticity. More- 
over, due to phase space restrictions we find that the 
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residual short-ranged interactions asymptotically decou- 
ple at low energies. The resulting RFL phase supports 
both gapless charge and spin excitations with no bro- 
ken spatial or internal symmetries, just as in a conven- 
tional Fermi liquid. But, due to the vortex sector of the 
theory, the RFL phase is demonstrably a non- Fermi liq- 
uid, with a gapless "Bose surface" of Rotons and with 
ODQLRO in the Cooper pair field but no Meissner ef- 
fect (see below). Moreover, the quasiparticles at the RFL 
Fermi surface are sharp (in the sense of the electron spec- 
tral function), but electrical currents are carried by the 
(quasi-)condensate. Even with impurities present the re- 
sistivity vanishes as a power law of temperature in the 
RFL. The power law exponent, 7, varies continuously, 
but is greater than or equal to one. For 7 < 1, the zero 
temperature RFL phase is unstable to a quantum con- 
finement transition, which presumably drives the system 
into a conventional Fermi liquid phase. 

To describe the RFL, we start with the general La- 
grangian, Eqs. H68I71|I . and first make the same approxi- 
mation as in the RL of expanding the roton hopping term 
to quadratic order. That is, to leading order, Ckin ~ ^r, 
with 

K 1 
= -^[dxyO - -;^{d^{ay - ay) + dy{a^ - a^)}f 

+ yN5.(a,-a,)]^ (172) 

As before for the RL, this approximation will be corrected 
perturbatively by "charge" and vortex hopping terms, 
which will not be expanded. 

Turning to the fermionic sector, we assume for the 
moment that the power-law decay (ODQLRO) of the 
Cooper pair field is sufficiently rapid that the pair field 
term Aj can be neglected. We argue later this is cor- 
rect for Ac > A* > 3/2. This gives a non-anomalous 
fermionic Lagrange density Lf, which we further presume 
is well-described by Fermi liquid theory, again checking 
the correctness of this assumption perturbatively in the 
couphngs to /3p and e^. Hence we replace £/ ~ £° -I- £j, 
with (working for simplicity at zero temperature) 

4 = /rV(5o~/i)/r.-tE/r+x,-/r., (173) 

3 

L) = -z/3o(r)/t./„ (174) 
j 

+ E[|/5.W' + |(/3,(r) +7re,(r) + A,(r))2]/rV^^.,/r,. 

3 

Note that, to leading order, the fermionic dispersion is 
controlled by the sum of the two hopping amplitudes, 
t = ts + te. Here we have included explicitly the physical 
external vector potential, Aj (r) , in the electron hopping 
term. Some care needs to be taken when treating Aj (r) . 
The above form is correct provided Aj (r) is also coupled 
into Cj in the quadratic Hamiltonian in the "canonical" 



fashion ncj — > ncj + Aj^ in the roton/plasmon portion of 
the Hamiltonian. It is not correct if this canonical cou- 
pling is removed by shifting e^, which is the procedure 
needed to generate the coupling in Eq. H74|l . If the 

latter form of the Lagrangian is used, the vector poten- 
tial should be removed from the electron hopping term. 
Either choice is correct if used consistently. 

The full Lagrangian that we then use to access the 
RFL phase is given by, 

Crfl = Ca + Co +J^r+^cs (175) 

with the definition, 

Co^^{doe~ao + aof. (176) 

The interaction terms between the Fermions and the 
fields Pj and Cj in will be treated in the ran- 
dom phase approximation. Doing so, one arrives at a 
tractable, if horrendously algebraically complicated La- 
grangian describing the RFL, which is quadratic in the 
fields 9,a^,ej,a^ and P^. This makes calculations of 
nearly any physical quantity possible in the RFL. 

Before turning to these properties, we verify (in the 
remainder of this section) the above claims that the cou- 
pling of fermions and vortices does not destabilize the 
RFL - i.e. it neither modifies the form of the low energy 
roton excitations at the "Bose surface" nor the fermionic 
quasiparticles at the Fermi surface. 



A. Quasiparticle Scattering by Rotons 

In this subsection, we show that the coupling of the 
electronic quasiparticles to the vortices does not destroy 
the Fermi surface. To do so, we will integrate out the vor- 
tex degrees of freedom to arrive at effective interactions 
amongst the quasiparticles. This procedure is somewhat 
gauge dependent. To provide a useful framework for the 
calculation of the electron spectral function in the follow- 
ing section, we will choose the gauge V ■ P — neijOiCj, in 
which the /, operators create fermionic quasiparticles 
with non- vanishing overlap with the bare electrons, with- 
out the need for any additional string operators. This is 
essentially equivalent to working in the electron formula- 
tion of Appendix ^ 

Since this gauge choice explicitly involves eijdiCj, we 
must employ a path integral representation in which the 
transverse component of the electric field, e*, has not 
been integrated out. It is further convenient to fix the 
two remaining gauge choices according to V • a = = 0, 
and to integrate out the field cq in the uq limit. 
The full RFL Lagrangian density (before imposing the 
constraint /?' — ire^j) then takes the form, 

vort H" 

{(3^) + £}{(3^), (177) 
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with 



and we have the definitions, 



where 



2u„ 



H — eijai{dj(3o - do(3j) + C0{aj - a*), 



+ !^[A,ye+^id,ay + dya,)]^ (179) 

To assess the perturbative efi^ects of the vortices upon 
the electronic quasiparticles, we wish to integrate out the 
vortices perturbatively in the coupHng of the gauge field 
to the fcrmions. For this we require the correlation 
functions of the (3^ fields (neglecting the couplings inside 
£/, and in particular to lowest order just {(3fj,(3iy)). To 
obtain the latter, we add a source term to the Lagrangian, 



+ i(Aoi/3o + A-/3). 



(180) 



Here we have included an extra factor of i with /3o to 
compensate for the factor of i present in the coupling 
of /3o to the fermion density. Upon fully integrating out 
the e*, 0, a*, ttj fields, the coefficient of A^Ai, in the 
effective action will give (half) the desired correlator. We 
perform the integration in two stages. First, imposing 
the constraint /?' — ire*, we eliminate e*, shift aj — *■ 
aj — {dj9 — a*j), and integrate out the 9, a*-, and fields. 
One obtains Cvort C,vort{otj, A^), with 

2 

Cvort = -^{etjdiaj + iTrXo)^ + -^[{d^ay)'^ + {dya^)^] 
+ ^^((9oai - Tre^jAj)^ - ty cos{aj), (181) 

with Ao = Ao — tp- 

In Eq. (|181(l we have added back in the vortex hop- 
ping term, — —tyCOs{aj), neglected in the RFL La- 
grangian. In the RFL phase, the vortex hopping is ir- 
relevant, and scales to zero at low energies. If we put 

— 0, the remaining integrations over aj are Gaussian 
and can be readily performed. We will return to the 
effects of non-zero vortex hopping upon the fermions in 
Sec. lVIlSI The final effective action then takes the form. 



S,^ff{X^) = -\ f C/^"J(k,c^„)A^(k,c^„)A.(-k,-a;„). 

(182) 



Here 



9 / A"- 



(183) 

specifies the propagator: (/3o/3o)o = C^o? , {Pil3j)o = 
~Uij\{iPof3j)o = —Uq^\ where the superscript zero re- 
minds us that this is the result to zeroth order in ti, = 0, 



uoo = + "-^Lol + (184) 

Uxx ~vl[vl\JCxK.y\^ +v\\JCx\^uj'i], (185) 
^vv = -vl[il\JCx]Cy\^ + vl\JCy\^u;l], (186) 



(187) 

viICxitOn[2io'i + vlilC^ - IC^)]. (188) 



2 r- . ro, ,2 



In the d.c. limit, these interactions simplify consider- 
ably, and one obtains the simple results 



C/^^^(k,C.„=0) 



4 \vo 

2 



;7^)(k,c^„ = o) = 

Uv 



(189) 
(190) 



andC/^°^(k,o^„=0) = 0. 

Since i/9o and /3 couple to the fermion density and 
current respectively, Seff mediates an effective fre- 
quency and wavevector dependent interaction between 
fermions. Since, in the d.c. limit, — ?7QQ''(k, w„ = 0) 
and C/j^°'' (k, LUn — 0) are finite and wavevector indepen- 
dent, they describe local (on-site) repulsive quasiparti- 
cle density-density and attractive current-current inter- 
actions, respectively. 

Generally, a repulsive density-density interaction be- 
tween fermions will lead to Fermi liquid corrections in 
the quasiparticle dynamics and thermodynamics, but will 
not destroy the Fermi surface. On the other hand, the 
current-current interaction for near-neighbor quasiparti- 
cle hopping can be rewritten in terms of antiferromag- 
netic, near-neighbor repulsive, and pairing interactions: 



H, 



-ni ni, 



2AfA{, 



(191) 



up to a shift of the fermion chemical potential, with S = 
P^f , (with a the vector of Pauli matrices), = f^f, 

Af — f^fi and — fffj- Here, J ~ t'^/uy is inversely 
proportional to the vortex core energy. One expects that 
the antiferromagnetic interaction can lead to an a Cooper 
instability in the d-wave (or extended s— wave) channel. 
The repulsive pairing interaction interaction dis-favors an 
s-wave Cooper instability. Hence it seems possible that 
for small Uy, for which this J is large, a spontaneous 
quasiparticle pairing instability may occur, most proba- 
bly of d-wave symmetry. Another possibility, which ap- 
pears very natural for extremely small doping x — > 0^ 
(and Uy — > O"*"), is that the antiferromagnetic interac- 
tion drives an antiferromagnetic spin-density-wave insta- 
bility. We will discuss both possibilities in the discussion 
section. 
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B. Roton Scattering by Quasiparticles 

Having established that the Fermi surface remains in- 
tact (apart from a possible BCS-type pairing instability) 
in the presence of gapless rotons, we need to see how the 
fermions feed back and effect the roton modes. To this 
end, we will integrate out all fields except 6, treating the 
fermions within the random phase approximation (RPA), 
to study the effects upon the roton dispersion. The RPA 
is complicated by the two distinct amplitudes, tg and ig, 
describing spinon and electron hopping processes, the lat- 
ter coupling to the electric field ej as well as the (ij gauge 
field. We begin with the RFL Lagrangian in Eq. (|175|l . 
It is convenient to first shift — > -I- then integrate 
out ao and take uq 0, which constrains V x /? = ttV • e 
and ao = —do6. Choosing in addition V • /3 = — ttV x e, 
we essentially return to the electron formulation, with 
f3i ~ — Trei. This choice is convenient, since the f3i and 
TTCi fields present in the electron hopping term precisely 
cancel, and only the — Tre^ appears "like a gauge field" in 
the spinon hopping term. The fermionic quasiparticles 
can then be integrated out in the random phase approx- 
imation (RPA). One finds 5*/ — > Srpa, with 

Srpa{I3^) - i / [noo|/?o(k,c^„)P + ^nn|eKk,c^„)P 

+^2i,(l-y)r„„|e(k,c^„)|2], (192) 

where T„„ = {f^^f^_^_^^) is the nearest neighbor kinetic 
energy, and t = tg+te- Here Hqo and Hn are respectively 
the density-density and current-current response kernels 
(polarization bubbles -I- diamagnetic contribution in the 
case of Hii) that would obtain for the Fermi sea were 
a single gauge field minimally coupled to the fermions. 
These depend upon the band structure. We will not 
require particular expressions for these quantities, but 
will use the fact that both Hoo(k, w„) and Hii(k, a;„) 
are finite and generally non-vanishing at Gxed wavevec- 
tor |k| > and tOn = 0. Since we will focus upon the 
low energy rotons, which occur near the principle axes in 
the Brillouin zone, we have dropped terms in Eq. H192|) 
that are small for ujn ^ <^ ky and tOn ^ ky <ti k^ 
(e.g. due to the parity symmetry, x —x, of the 2d 
square lattice, Hoi(/cx — 0,ky,LUn) vanishes). We also 
note that Eq. (|192|l does not have the usual RPA form 
even in this limit, due to the non-minimal coupling form 
of the fermion Lagrangian (minimal coupling is restored 
only for ts —^t). 

It may be helpful to keep in mind the forms for a circu- 
lar Fermi surface (e.g. valid at small electron densities), 
where, at small wavevectors and frequencies, one has 

Hoo(k,c^„)~m*, (193) 

where m* ^ is the effective mass and 

Hn(k,c.„)~^[k2 + i4], (194) 
m* vpk 



which is valid for |w| < vpk. We stress, however, that 
our results do not depend upon these particular forms. 

Since Srpa is quadratic, one can perform the remain- 
ing integrals straightforwardly. We choose V • a ~ V • o? — 
0, and integrate out a, /3o and e, to obtain 

S = [ f ^2 I K,r ^vl+Vl ^2 , ^n^^ g,2 

+ f |/C./C,^?+i^^^ap}, (195) 

where u„ = Uy+n'^ts{l~ts/t)Tnn,Uy = u^-t-vr^t^/t^Hn, 
and vl = v\ + KrTT^HoUQ. Without loss of generality, we 
focus upon the branch of rotons with ujn ^ k^ ^ ky ^ 
0(1), for which the first term in Eq. (|195ll is negligible, 
and the remaining a integral can be carried out to obtain 
finally the effective action 

Srot = \l P H + vlot{ky)kl] \0\\ (196) 

with 

We have thereby shown that even a gapless Fermi sea 
does not lead to a qualitative change in the gapless Bose- 
surface of rotons. Due to the fc^-dependence of Hqo and 
Hii (implicit in vi and Mi,), the roton velocity is now seen 
to depend upon ky. Additional "direct" quasiparticle- 
quasiparticle interactions (not mediated by the rotons) 
would in any case similarly renormalize Vrot ■ But the lo- 
cation of the Bose surface and the qualitative low energy 
dispersion of the roton modes are seen to be unaffected by 
the fermions. Together with the earlier demonstration of 
the stability of the Fermi surface, this result establishes 
the RFL phase as a stable 2d non-Fermi liquid with both 
gapless charge and spin excitations. 

VI. INSTABILITIES OF THE ROTON FERMI 
LIQUID 

As for the RL, the RFL has potential instabilities to 
superconducting and Fermi Liquid states, driven by ef- 
fects/terms neglected in the previous subsection. We ad- 
dress each in turn now. 



A. Landau Fermi Liquid Instability 

The arguments of Sec. IIV Al for determining the rele- 
vance or irrelevance of the neglected vortex hopping term 
within the simpler RL continue to hold for the full RFL, 
provided the proper renormalized roton liquid parame- 
ters are employed. In particular, the vortex hopping term 
continues to be described by a 1 -I- 1-dimensional scaling 
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dimension A^,, which is, however, renormaHzed by the 
statistical interactions with the fermions, to wit: 



1 

= — 
An 



27r Vrot{ky)' 



(198) 



with Vrot{ky) given from Eq. (|197|l . The same arguments 
thus continue to apply, and the vortex hopping term is ir- 
relevant for Ai, > 2. For A„ < 2, we expect an instability 
to a state with proliferated vortices. The dual Meissner 
effect for this vortex condensate confines particles with 
non-zero gauge charge, as discussed in Sec. IIV Al Com- 
ing from the RFL, the natural expectation is then that 
the system becomes a Fermi liquid. This hypothesis is 
fleshed out in Appendix IXI 



B. Superconducting Instabilities 

As for the RL, the RFL can also be unstable to a 
superconducting state. However, the inclusion of the 
fermionic quasiparticles opens new routes to supercon- 
ductivity from the RFL. We explore each of these in turn. 



1. "Charge" hopping 

As for the vortex hopping considered above, the ar- 
guments for the relevance of the "charge" hopping in 
Sec. IIV B 21 for the RL continue to apply with only a 
renormalization of the roton liquid parameters. In par- 
ticular, the scaling dimension defined from the charge 
hopping, Eq. (|169|l . is modified to 
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With this modification, the charge hopping processes be- 
come relevant for Ac < 2, as before. It should be noted 
that, including the renormalizations due to scattering by 
fermions, AcA^ ^ 1/2, owing to the ky dependence of 
the Vrot and Uy . 

As we established in Sec. HVB2l bv employing the pla- 
quette duality transformation, when the charge hopping 
processes are relevant the rotons become gapped out and 
the system exhibits a Meissner effect. Here, we briefly 
comment on the corresponding fate of the fermionic 
quasiparticles, which are gapless across the Fermi surface 
in the RFL phase. The relevance of the charge hopping 
in the plaquette dualized representation, indicates that 
it is not legitimate to expand the cosine term that enters 
into the roton hopping Lagrangian, Cr- Rather, in the 
original vortex representation, the state corresponds to 
a "vortex vacuum" , and the properties of the state can 
be accessed by taking all of the vortex hopping processes 
small, tv,t2v,Kr <^ Uy. At zeroth order in the vortex 
kinetic terms, the full Hamiltonian of the U{1) vortcx- 
spinon field theory in Eq. I|43|l is independent of aj, so 
that f3j commutes with H and naively can be taken as 



c-numbers. A simple choice consistent with the condi- 
tion TVr = is Pj — 0, and in the vortex vacuum one also 
has that ej = 0. In this case, the full fermionic Hamilto- 
nian describing the quasiparticles in Eq. I|48|) reduces to 
the Bogoliubov-deGennes form, apart from the coupling 
to e* which describes the Doppler shift couplings to the 
super-flow. In particular, the string operator which en- 
ters into the explicit pairing term equals unity, 5^ = 1, so 
that the fermionic quasiparticles are paired. We denote 
by |0) / the ground state of Hj with f3j = Cj = 0, which 
is easily found by filling the Bogoliubov-deGennes levels 
below the Fermi energy. Our naive ground state is then 



>o = |0)/®|/3, =0)®|A^, = 0). 



(200) 



In the U(l) sector, unfortunately, we must take some- 
what more care, since the condition (3j = is in fact 
inconsistent with the gauge constraint Qf{Ar) = 1. We 
can, however, project into the U(l) subspace using the 
operator 

^ = n ^ifrJra - -^r^^jir - w)). (201) 
r 

Since [P,H] = [P,^i,(Xr)] = 0, the state 

|0) = F|0)o (202) 

satisfies all gauge constraints and remains an eigenstate 
of H (with zero kinetic terms). This establishes that the 
resulting superconducting state is an ordinary BCS-type 
superconductor with paired electrons. 



2. BCS instability 

As we saw in Section V, the coupling of the spinous to 
the gapless rotons in the RFL phase leads to Heisenberg 
exchange and pair field interactions between the fermions 
with strength J ^ t1/uy. In the physically interesting 
limit t « is 3> te, the quasiparticles move primarily as 
spinous, i.e. without any associated electrical charge, and 
hence do not experience a long-ranged Coulomb repul- 
sion. Thus one may imagine that the above interactions 
could lead to a BCS pairing instability at "high" energies 
(still below the quasiparticle Fermi energy to make the 
Fermi liquid description appropriate, but quantitatively 
large compared to e.g. more conventional BCS critical 
temperatures) of order J. Here we explore the properties 
of the resulting phase which emerges when the fermions 
pair spontaneously and condense. To this end, we focus 
on the fluctuations of the phase, $ of the fermion pair 
field, ififi) = A/e**. Keeping A/ fixed, and working 
once more in the electron gauge Pi — — Trci, we integrate 
out the fermions entering in Cj^pl to generate an effec- 
tive action for ^,/3oi and e^. Specifically, one obtains 
C°j- £} -> £$, with 



.9/ 



(9o$ - 2/3o) 



(203) 
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Here gf is of order the density of states at the Fermi 
surface, Ai is (time-independent for simphcity) exter- 
nal vector potential, and p( is the fcrmion pair ("su- 
perfluid" ) density. The final term represents the low but 
non-zero temperature corrections to the superfluid den- 
sity appropriate to the case of d-wave pairing. Here a 
is a non-universal constant of order one representing the 
effects of Fermi liquid corrections, or equivalently, high- 
energy renormalizations of the "Doppler shift" coupling 
constant. 

The excitations and response functions of the system 
can then be obtained from the RFL Lagrangian by shift- 
ing a^ + a^ and then integrating out ao, aj and uq. 
Having made the replacement, Cf one thereby ob- 
tains Crfl — >■ ^eff, with an effective Lagrangian given 
by Cef f — Crot + Cq, , where 
u 1 

I TT 



(204) 



+ '^{eijd^ajf + '^[A.^yO + i(9^ay + dya,j,)f . 

First, it is instructive to consider the transverse elec- 
tromagnetic response, in particular consider a static, 
transverse external gauge field, d^Ai =V-A — 0. Since 
the external gauge field is at zero frequency, aj,(3Q, 9, 
and et are decoupled from Ai and ei in this limit, and 
we may thus neglect all but the first term in Eq. 12041 In 
addition, <i> is decoupled from e; (et) as well, so we may 
simplify the effective Lagrangian to 



^eff = ^(e/ + — )" + 2tep{A2 + 2n%plef 

I TT 



(2^i,e/-2teAt)2. (205) 



Integrating over e; then gives the superfluid stiffness 
Ks as the coefficient of 2A^ (corresponding to the pair- 
field phase stiffness) in the effective action: 

K, = (206) 
p{{tu^ -f iirHetspO (tuy + AirH^tspif 2akBT 



+ 4:TT^tspl 



(u, + AirHspiy tAf 



to linear order in temperature. Note that, for p( ^ 0, the 
system is a true charge superfluid, displaying a Meissner 
effect. 

We would next like to establish the fate of the roton 
excitations. To do these, we let Ai = 0, specialize to 
T = 0, and consider the appropriate limit w„ ~ fc^, ^ 
ky ^ 0(1), integrating out fields to obtain an effective 
action for B alone. Integrating out $, one obtains in this 
limit 



using the above conditions on u;„ and k, and doPo <C 
eijdiCj, which follows in this limit. Further taking the 
gauge V • (5 — 0, one sees that ej couples only to doa 
in Eq. (|204(l . and thus generates only negligible terms 
at low frequencies and may be dropped. We therefore 
need only the effective action for a ~ at, 9, (3o and e;, 
^eff — ^eff, which in this limit becomes 

«e// - 7. 6/ + ^-5-2 Po - |^-y|w„e/y 

^|M4-^aP(208) 



+ 



2 



Performing the integration over a, [3q and e/, one obtains 
the final effective action for 9: 



effiO) 



1 

2 U,: 



^TT^pit, 



{^l + v^tkl) \0\M209) 



where 



„--.2 



u,{vl + ^4) + 2^^gfvlvl 



(210) 



Thus, despite the superconductivity induced by quasipar- 
ticle pairing, the gapless roton excitations survive, with 
some quantitative modifications to their velocity and cor- 
relations. 

As we shall explore further in the concluding sections, 
in the limit of a very small "bare" vortex core energy 
appropriate in the under-doped limit, ~ a; — > 0, there 
is a large energy scale for pairing, J ~ tj./u^. It is then 
natural to assume tpf. m t{l — x) ^ Uv If we presume 
that the fermionic kinetic energy is predominantly due 
to "spinon hopping" , t tg ^ t^, then one has 



^ + t,pf-i^+UpfrJ^^ (211) 



(207) 



Thus, despite the large bare superfluid stiffness coming 
from the BCS pairing of the quasiparticles, the renormal- 
ized stiffness is small, set by the bare vortex core energy, 
Ks — ^ x or the electron hopping (presumed 
small). It is the renormalized stiffness which determines 
the vortex core energy, and sets the scale for the finite 
temperature Kosterlitz-Thouless superconducting transi- 
tion, Tkt ^ Uy ^ X. In this way, one can understand the 
large separation in energy scales between the pseudo-gap 
line at scale J and the superconducting transition tem- 
perature, Tkt ~ X. Unfortunately (since it is in appar- 
ent conflict with the small amount of experimental data 
for this quantity), along with this small superfluid stiff- 
ness, one obtains a small linear temperature derivative, 
dKs/dT\rp^Q, due to the same mechanism. This is sim- 
ilar to results for the U(l) gauge theory for the t — J 
model. 

To complete the analysis, one should consider the ef- 
fects of the heretofore neglected explicit pairing term 
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in this novel "rotonic" superconductor. In particular, 
one may imagine that, once the fermionic pair field 
has condensed, it may induce true off-diagonal long- 
range order in the rotonic sector through the proxim- 
ity effect. Naively, ODLRO appears incompatible with 
gapless rotons, so one may expect the explicit pairing 
term to induce a gap in the roton spectrum. While 
this is possible, it is easy to see that it is not in- 
evitable. To see this, note that, in the rotonic super- 
conductor, the fermionic pair-field in the explicit pair- 
ing term may be replaced simply by its mean-field value, 

AiCr+xio-eo-o-'Crcr' = l^lCr+^^cri-aa' C-ra' = A/. This leads, 

in the roton sector, to an additional term, 

i^A = -A/^(<r+x, +h.C.)- (212) 

Indeed, Eq. H212() . since it embodies a linear coupling to 
the bosonic pair field, i3r,r-i-xj , will induce ODLRO in the 
B-c operators. However, this need not itself be a mecha- 
nism to induce a roton gap. Noting that from Sec. lIII CI 
the boson pair field correlators are power-law in form, we 
expect that for Ac sufficiently large, perturbation theory 
in A J will be regular and convergent, and the gapless 
rotons will be preserved. Due to the anisotropic struc- 
ture of these correlators, we cannot reliably determine the 
relevance or irrelevance of Aj by simple power-counting 
arguments. We note that for Ac > 3/2, the induced 
ODLRO is expected to be linear in A/, i.e. the bosonic 
pair-field susceptibility is finite. However, the criterion 
for irrelevance is probably more stringent, e.g. Ac > 3. 
To determine the precise value A* such that for Ac > A* , 
the explicit pairing term remains irrelevant after conden- 
sation of quasiparticle pairs, would require a more careful 
analysis of the structure of the perturbation theory in A j . 
We leave this for braver souls, and content ourselves with 
the observation that there is a range of stability to this 
perturbation. We note that, however, since the RL itself 
even without the pairing term is always unstable to either 
vortex or "charge" hopping, the rotonic superconductor, 
with true gapless rotons, exists at best as an interme- 
diate energy scale crossover. Nevertheless, provided the 
inevitable roton gap is small (i.e. for small A/, tc, we 
expect the gap onset will only slightly modify the val- 
ues for the superfluid stiffness and its linear temperature 
derivative determined above. 



3. Explicit pairing term 

The final potential instability of the RFL we consider 
is from the explicit (spinon) pairing term A. Recall from 
Sec HTTTTl 

Ha - ^e*'3^('^U,[5r]Vr+x,.W/ra'+h.C.(213) 

3 

= ^ AjB^_^^^^.Cr+x,,Te^CT'Cro-' +ti-c., (214) 
3 



where the last line is written in electron variables, or 
equivalently in the "electron gauge" with j3j = —TfCj. 
This representation is convenient because it eliminates 
all unphysical gauge fiuctuations, which, although they 
do not appear in any physical properties, may enter inad- 
vertently through approximations. Had we considered in- 
stead a Hamiltonian with s-wave pairing, we would have 
had 

^ A^[5r]Vr.e..'/r.'+ll.C. 

A"Stcr^e„^,Cra, -^h.c. (215) 

Since we have already determined the correlations of the 
boson pair field operator Br (which exhibits ODQLRO), 
in Sec. IIII CI all the operators appearing in Eq. H215|) 
have well understood properties at this stage, and so we 
will discuss this case for simplicity. For the physically 
more interesting scenario of d-wave pairing, we require 
instead the behavior of the bond pair field correlations. 
As discussed in Sec. lIII (T3l this behavior is qualitatively 
identical to that of the local pair field, with a possible 
renormalization of Ac. Hence we believe that the results 
we obtain for local (s-wave) pairing - in particular that 
the pairing term is irrelevant for sufficiently large Ac - 
carry over to the c?-wave case, provided a possible (and for 
the moment unknown) 0(1) modification of Ac is made 
in the relations. 

We would like to determine the "relevance" of TYa in 
the renormalization group (RG) sense, i.e. whether its 
presence destabilizes the low energy properties of the 
RFL. Unfortunately, due to the extremely anisotropic 
nature of the rotonic spectrum, and the very different 
nature of the low energy electronic quasiparticle states 
at the Fermi surface, we do not know how to formu- 
late a proper RG transformation. However, we do note 
that correlations of TYa decay as power laws in space and 
imaginary time in the theory with A = 0, since then these 
correlations factor into G"^p (with power-law behavior de- 
scribed in Sec. lIII G|) and the fermion pair-field correlator, 
which has the power law form characteristic of a Fermi 
liquid. Clearly, for sufficiently large Ac, the Cooper pair 
propagator G"^P(r, r) will decay sufficiently rapidly that 
perturbation theory in A does not generate any (primi- 
tive) singularities. However, determining the critical A* 
above which this occurs is beyond the scope of the cur- 
rent study. A simple argument clearly bounds A* > 3/2. 
In particular, one may attempt to integrate out the vor- 
tices perturbatively in TYa in a cumulant expansion. The 
first non-trivial term in this expansion occurs at second 
order, and generates an attractive fermion pair-field to 
pair-field interaction whose vertex is simply the Fourier 
transform of G"^P(r, r). For Ac < 3/2, a simple scaling 
analysis indicates that this Fourier transform is divergent 
at q = = 0. Thus for such values of Ac, this attrac- 
tive Cooper channel interaction will overwhelm any other 
repulsive interaction that might be present at low ener- 
gies, leading to a Cooper instability and pairing. This 
analysis, however, neglects higher cumulant terms which 
are certainly present in integrating out the vortices, and 
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which are presumably crucial in determining the ultimate 
limits of stability of the RFL. Nevertheless, and this is all 
we shall require at present, it is clear that A* exists and 
is not infinite, so that a non-vanishing region of stabil- 
ity also exists. When Ac < A* the explicit pairing term 
will be "relevant" , destabilizing the RFL phase, presum- 
ably driving it into a conventional superconducting phase 
with paired electrons and gapped rotons. Given our lack 
of knowledge of A* , the additional uncertainty in the de- 
cay exponent of ODQLRO is not particularly damaging. 



VII. PROPERTIES OF THE RFL PHASE 

Having established the existence of the RFL phase, we 
now discuss some of it's properties. Since the rotons are 
effectively decoupled from the fermions at low energies, 
much of the physics of the RFL phase follows directly 
from the results we established for the charge sector in 
Sections III and IV. Specifically, one expects three dif- 
ferent gapless excitations in the RFL phase - a gapless 
longitudinal plasmon, a Bose surface of gapless rotons 
and a set of particle-hole excitations across the Fermi 
surface. The spin physics in the RFL phase is then qual- 
itatively similar to that of a conventional normal metal. 
Also like a metal, the RFL has ODQLRO in the Cooper 
pair field. Unlike an ordinary metal, however, this power 
law off-diagonal order exists with two distinct unrelated 
powers, one originating from the two-particle excitations 
of the Fermi sea, and the other from the gapless rotons. 
Because the rotons and electronic quasiparticles exist as 
nearly independent excitations in the RFL, one would ex- 
pect a sharp electron spectral function despite the critical 
rotons. While this is true, as we demonstrate below in 
detail, the quasiparticles do scatter appreciably and in 
a singular manner from the gapless rotons at particu- 
lar "hot spots" on the Fermi surface. Nevertheless, even 
at these hot spots, the decay rate vanishes more rapidly 
than the electron's energy. 

Despite the existence of such long lived electron-like 
quasiparticle excitations, the electrical transport in the 
RFL phase is strikingly non-Drude like, as we now 
demonstrate. Specifically, in the presence of impurities 
the electrical conductivity at low temperatures is domi- 
nated by the quasi-condensate in the charge sector, di- 
verging at low temperatures as cr '--^ J'-t -^ith 7 > 1. 
But as we shall see, the Hall conductivity behaves very 
differently, being dominated by the electron-like quasi- 
particles. 



A. Electrical Conductivity 

Here we evaluate the electrical resistance in the RFL 
phase. As in Section IVA, it will be necessary to include 
the effects of the vortex hopping term, even though vor- 
tex hopping is technically "irrelevant" when 7 > 1. The 
RFL phase exhibits the same "emergent symmetry" as 



in the Roton liquid - the number of vortices on every 
row and column being independently conserved - so in- 
clusion of vortex hopping is necessary to generate any 
dissipative resistance whatsoever. In order to access the 
Hall conductivity we apply a uniform external magnetic 
field. We choose a gauge with ^0 = 0, and set. 



A,{r,r)=Af{r)+A,{r), 



(216) 



with tijdiA^ — B the external magnetic field and Aj a 
time-dependent source term used to extract the conduc- 
tivity tensor. As might be expected it is important to 
include the effects of elastic scattering from impurities. 

It is convenient to employ, as in Sec. IVI B 21 in which 
the fermions have been integrated out, their effects felt 
only through an additional contribution to the effective 
action. In particular, we choose, as in Sec. IVI B 21 the 
electron gauge, and write 



If 



trot - tv cos{di9 + a-j) -I- Sf{e,,Po,A,), 



(217) 

with Lrot given in Eq. (|2U4|I . Here Sf represents the 
fermionic terms in the action. We proceed by shift- 
ing Cj — > Bj — A^ /tt, which takes Aj Aj in Lrot 
without any further changes since A^ is time indepen- 
dent. This has the effect for the fermions of making 
the magnetic field appear uniformly in both electron 
and spinon hopping terms. Integrating out the fermions 
then pcrturbatively in e^, /3o, and Ai effectively replaces 



S 



f 



where 



Sf = I^|;3o(k,c.„)P + ^T„„|e,f + |r„„|i,.p 

+ ^^{T^tsei + teAi)k^^^{TTtsej -I- ie4j)-(kL^„)- (218) 

In Eq. H218|l . both Hoo and Hy are functions of k and 
ujn (and B through A^ which appears in the electron 
hopping term), and we have neglected a cross-term Hoi 
between Pq and the spatial gauge fields, which is neg- 
ligible in all the limits of interest. We will require the 
behavior of H^ in two regimes. Since the external fields 
are spatially uniform, we will need H(k = 0,cij„) at low 
frequencies 



ny(k: 



0, w„) 



\uJn\Si 



(219) 



where a(j is the conductivity tensor for the fermions. Ro- 
ton fiuctuations are dominated in contrast by w„ ~ fc^, ^ 



ky ^ 0{1) (or the same with 



In this limit, H; 



becomes a non-trivial function of the 0(1) component 
of the wavevector, but has the useful property (due to 
square reflection symmetry) of decoupling in the longitu- 
dinal and transverse basis, 



Ilij(k,uj„) « Ht 



Sij - 



/C2 



+ Hz^. (220) 
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In this limit, we note that lit is related to XIn of Sec. IV Bl 

by IIii =flt+ Tnn- 

As discussed above, the conductivity is finite only once 
the vortex hopping is included to break the row/column 
symmetries of the RFL. Hence to extract the conduc- 
tivity, we must compute the effective action for Aj to 
O(i^), which gives the first non-trivial correction. This 
requires only Gaussian integrals, with no further approx- 
imations. However, for ease of presentation, it is conve- 
nient, analogously to Sec. IIV iTTl to take the simplifying 
limit Uy/vQ -\- Hgo <SC I/k^- In this limit the fluctua- 
tions of (3q are extremely strong, which in turn strongly 



suppresses the fluctuations of a*(k, Ci;„) except at k = 0. 
Furthermore, choosing the gauge V • d? = 0, we may 
thereby take Q;i(r, r) to be a function of r only. Do- 
ing so, we may drop all spatial derivatives of aj in Crot- 
Furthermore, since fluctuations of Uj are only temporal, 
it can be accurately treated in an RPA fashion. 

To carry out the RPA calculation, we first perform the 
integral over e^, which gives an effective action in terms 
of the remaining 9, ai, and Ai fields, Seff S'^^j = 
S[{A{) + S'^{a,,A,) + S'^{e,ai), with 



S[ = 
S' = 



I l^iilts + te)r„„Jy + 75 (ie + V^tl + 277teis) Hy (k = 0, W„ ) } A: ("Wn) , 



2 r 



1 



r]UJ. 
u„ 



t2 

^H,,(k = 0,L^„) 



ai{-UJn)aj{Un) - ''^^eijai{-UJn)Aj{uJn) 



+ a^), 



(221) 
(222) 
(223) 



where 77 = u^/ {uy -t- n'^tsTnn), and Uy = Uy + Tr^tgTn 



-KHlHt/t^ = Uy + T:H^,Uii/t^ as in Sec. |VB1 In 
Eqs. 1)2211) . (|222|l . L'^ is the system volume (number of 
sites in the square lattice), arising since ai,Ai are spa- 
tially constant. In obtaining Eqs. H221^ I222|l . we ex- 
panded to linear order in H^j (k = 0, cjn) since we are in- 
terested ultimately in the low- frequency limit, and Hy is 
linear in frequency. In Eq. H223I) . we used the fact that 9 
coupled to only to the longitudinal part of and hence 
only to Ht . Note that the quadratic part of Eq. 1)223(1 re- 
produces precisely Eqs. (|196l - ll97Tl . in the hmit vq ^ vi, 
as assumed herein. 

We now can integrate out 6 to 0{tl). This proceeds 
identically as in Sec. IIV A II with aj replaced by —aj, and 
with the renormalized roton liquid parameters Uy ^Uy. 
Hence we obtain the correction S'-^{9, ai) S'.^{aj), with 



C// 

•^3 



4 



'^(w„)|aj(a;„)|' 



(224) 



with 7^(w„) |t^„r+'' as in Eq. piT|) of Sec. IIV All 

but 7 — 2lS.y — 3 with the renormalized A^, given in 
Eq. fTM . 

With this replacement, the remaining quadratic inte- 
gral over ai can be easily performed to determine the 
physical response kernel, S'i + S2 + S'^ Sresp{Aj), with 



il,(-c^„)H«^^(zc^„)ij(L^„),(225) 



where 



H,7^(zc.„) « nf.+nr;*, (226) 



and 



H{ 



H''°* 



i2 



Uy ultlTZ{uJn) 



Hy(k = 0,w„),(227) 



(228) 



•2 2^2^2 

The conductivity is obtained from the Kubo formula as 



(229) 



As in Sec. IIV ATI wc see that the rotonic contribution to 
the conductivity is not perturbative (at low frequencies) 
in ty, so to capture the expected behavior, we replace it 

by, 



yrrot 
ij 



This gives 



(230) 



a.,(c.,T)=a[f(c.,T) + 



t2 



a[,(a;,r). (231) 



Notice that from Eq. (|231|l , the conductivity is the sum 
of separate fermion and roton contributions, and that 
the only effects of the fermions upon the rotonic piece 
is to modify the exponent 7 (or A^,) implicit in TZ{uj). 
Moreover, the fermionic contribution vanishes for t,, — 0, 
as expected on physical grounds since the spinon hopping 
term does not transport charge. 
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Consider first the dissipative diagonal d.c. (sheet) re- 
sistance, R{T) = Re[a~^], which at low temperatures is 
completely dominated by the rotonic contribution, given 
as in Eq. (|150|l . 



i?(T) = 



ttT 

Vrot 



(232) 



except with Vrot and the exponent 7 renormalized by 
the interactions with the fermions. We thereby arrive at 
the important conclusion that the resistance in the RFL 
phase varies as R{T) ^ , with 7 > 1. When 7 < 1 the 
RFL phase is unstable to confinement into a Fermi liquid 
phase, and right at the quantum confinement transition 
the resistance is linear in temperature. 

The Hall conductivity and hence Hall angle on the 
other hand, are dominated by the fermionic quasi- 
particles. Specifically, a^y ~ cr^j,, and in the d.c. 
limit the fermionic Hall conductivity can be approx- 



imately obtained from a Drude expression, 
ujcTf{nfe^Tf)lm ^ ^ct'j, where ojc = eB/ni is the cy- 
clotron frequency. The fermionic scattering time Tf has 
a variety of contributions, from elastic impurity scatter- 
ing to interactions with rotons, considered in the next 
section, and hence may (or may not) have temperature 
dependence of its own - in contrast to the diagonal con- 
ductivity, a XX ~ , which diverges as T ^ due to 
the rotonic contribution. These considerations suggest 
that the cotangent of the Hall angle, defined in terms of 
the resistivity tensor pij as cot(0//) = Pxxl Pxy^ will vary 
as. 



coiiQn) ^ — ~ 2 



1 
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The complete absence of a roton contribution to the 
Hall conductivity cr^^'" is a consequence of the magnetic 
field independence of the roton liquid Lagrangian, Cj^l, 
either in it's original or dual forms, Eqs. H9U|I and (|159|l . 
respectively. But as discussed in Sec. IVB, the dual rep- 
resentation of Crl allows for additional terms involving 
cosines of the dual field (/>, which are present due to the 
underlying discreteness of the vortex number operator. 
While being irrelevant in the Roton liquid phase, these 
neglected terms do depend on the external magnetic field 
and if retained will lead to a non-vanishing roton con- 
tribution to the Hall conductivity. But this contribution 
is expected to vanish rapidly at low temperatures. If on 
the other hand, these terms are relevant and drive a su- 
perconducting instability at low temperature, they will 
likely contribute signifigantly to the Hall conductivity. A 
careful analysis of the Hall response above the supercon- 
ducting transition temperature in this situation will be 
left for future work. 



B. Electron Lifetime 

Although the electron spectral function is expected to 
be sharp at zero temperature right on the Fermi sur- 



face in the ideal RFL "fixed point" theory, at finite ener- 
gies (or temperatures) one expects the electrons to scat- 
ter off the rotons through interactions neglected in the 
RFL, causing a decay. Here, we consider two contribu- 
tions to this electron decay rate at lowest order in the 
perturbations to the RFL. We consider two scattering 
mechanisms. First, scattering due to coupling of the 
quasiparticle current to boson currents through the /3p 
and Bi terms present in the fermionic hamiltonian. Be- 
cause singular bosonic current fluctuations are primarily 
induced by vortex hopping, non-trivial contributions to 
the fermion lifetime through this mechanism occur first 
at 0{t1). Second, we consider scattering of quasiparti- 
cles due to "superconducting fluctuations" , i.e. fermion 
decay mediated by the "Jospephson coupling" A to the 
bosonic pair field. Since we expect A ^ ^i, on physi- 
cal grounds (see Sec. IVHI B 2|) . this O(A^) contribution 
is naively much smaller than the former one. However, 
depending upon the parameters of the RFL, this need 
not be the case. This point will be returned to in the 
discussion. 

As discussed above, to compute the spectral function 
it is simplest to work in the gauge = 7re* , in which we 
may assume the electron operator Cro- ^ fra, and hence 
study 



(234) 



In other gauges, it would be necessary to include string 
factors as indicated in Eq. (|55|l . 

Neglecting the fluctuating /3p fields, one has 
Qo(k,u!n) — {i^n — Ck)""'^, with dispersion, = 
—2t cos(fcj ) — ps- Both fluctuations of the /3^i gauge fields, 
which mediate retarded interactions amongst the quasi- 
particles, as well as the explicit interactions in the pairing 
term, will induce a self energy I](k, w„), defined by 



5o(k,a;„) 



Ek - S(k, w„) 



(235) 



We will compute the lowest order corrections to the 
electron self energy, obtaining a single-particle lifetime 
from the imaginary part of its retarded continuation, 
Eret(k, (jj) = S(k,a;,i —^ —ilo -\- 0"*") in the usual way. At 
the level of this discussion, the two types of interactions 
act in parallel to scatter quasiparticles, so 



E(k,c^) = I],/(k,c^)-HSs/(k, 



(236) 



We will focus upon the imaginary part I]"(k, w), which 
describes broadening of the electron spectral function, 
l/T/(r) = S"(k,0;r). We have 



(237) 



In principle this single-particle "lifetime" is distinct from 
the momentum scattering rate which is relevant in dis- 
cussions of transport quantities. However, we will use 
the behavior obtained for 1/r/ as a crude guide to the 
quasiparticle transport as well, leaving a more careful 
treatment for future study. 
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1. Scattering mediated by vortex hopping- enhanced current 
fluctuations 

Taking into account quadratic fluctuations of the (3^ 
fields, the leading self energy correction due to current 
fluctuations is 



ImEret is given by, 



KfiKi^) = (241) 



E 



1 



VFfiVF,^—, 



with the definitions, VFj = dej^/dkj, vfo = 1- Here we 
have written the self-energy in terms of the full /3— gauge- 
field correlator, Ufj,i, to all orders in ty, rather than re- 
stricting to fi, = as we did in obtaining Eq. p83|l . 
Introducing a spectral representation. 



The self-energy can now be evaluated by considering 
, progressive orders in t^. To zeroth order, the expres- 

-C/,,.(q, LuJ, (238)g.^^g jj(o) ^^^^^ g^g^ JTHaiHi. 

We note that, because they come from the quadratic 
RL Lagrangian, they contain only simple poles. Fur- 
thermore, they are explicitly real functions of ia;„, so 
that, upon analytic continuation, their retarded correla- 
tor has zero imaginary part. Since ujj^J" = 0, one thus 
0. 



duj Ul^^{q,uj) 



TT U! - 



(239) 



with U';,{u) = ImUp{uj) = ImU^,{uJn ^ -zw + 0+), 
allows one to analytically continue to obtain the (re- 
tarded) self energy. 



has E^y"(k,w) 

Thus there is no broadening of the quasiparticle peak 
at zeroth order in the vortex hopping. The first non- 
trivial contribution to the quasiparticle lifetime occurs 
though the 0{t^) corrections to the gauge field propaga- 



tor Uuu 



This correction is obtained 



I],et(k,w) = S(k,c^„)|,^„^^+,o+- (240) 
For positive frequencies, (jJ > 0, the imaginary part, E" — 



by integrating out aj to second order in t„ starting from 
Eq. p81|l to obtain Seff{Xf_i) to O(t^) using the cumu- 
lant expansion. We shift ai(r, t) Q!i(r, r) — Vi{r,T), to 
eliminate the linear terms in ai in the Lagrangian, with 
Vi linear in given explicitly by. 



Va:i'k,UJn) = ' 2\f 2 , 2~T 1^0^^^!/ (""o^^ •^O + «Wn A^; ) - (w^ + W+ j/C^: P) (Wq Aq -f iw„ Aj, )] , (242) 



TT 



Wy(k,W„) = — ^— 2\/ 2 , 2~^[Vo)Ca:ICyivQlCyXQ + iUnXy) - {UJ^ + V^\JCy\ ){vQ)Ca:Xo + iiLlnXx)]- 

After this shift the correction to the effective action can be formally written, 

^S/ = -|EE/ ^{cos{a,{r,T) -v,{r,T))cos{a,{r,T) - v,{r,T)))^, (243) 

ij r.r' 

where the subscript a indicates a Gaussian average over with respect to the quadratic terms in Eq. H181|) . Since 
Vi is linear in A^, this correction is not quadratic in A^, indicating that the /J^ fiuctuations are not Gaussian. To 
evaluate the leading order self-energy correction, however, we require only the two-point function of (3^, and hence 
may expand to quadratic order in u^. One finds 

Si)\ = J (^\Vxik,U„)\^n{ky,iUr,) + \Vyik,Un)\^nik,,LUn)y (244) 

Here TZ{k,ujn) is the two-point function of the vortex hopping operator considered in Sec. lIV.^ i.e. 

n{k,ujn) = /(l-cos(fcj; + cj„T)){Je"«(^'^'^)e-*"''('^'^'°)^ 



^\'-'v) ( 2 I 2 l,2\A„-l 
l7r(A„ — ij 



sm-; 

r 



(245) 



where the latter behavior gives the leading non-analytic term for small ujmk, with A{/S.^) a constant prefactor. 
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An analytic quadratic term is also present (and larger 
than the given form for A„ > 2), but does not contribute 
to the imaginary part of the self-energy for the same rea- 
sons described above for the = terms. 

Next, we must analytically continue to obtain the 
U^Ji^"{'k,Lu). From Eqs. one immediately sees 

that u['^*(k, is purely real. Hence, the imaginary 
part comes entirely from the analytic continuation of 
R{ki,ujn). One has 



f)^^-^&{\uj\~Vrot\h\)sgn{uj). 

(246) 

From inspection of Eq. (|244|l . one thereby sees the 
imaginary part of the retarded gauge field correlator, 

Ujj,J (k, uj) , is the sum of two terms, which are non- 
zero only for |a;| > i^rot or w > Vrot\ky\, respectively. 
Thus the momentum integrals in the expression for E"^ 
in Eq. (|241|) . are constrained not only by < ek+q < ^ 
but also either ek+q < u}-Vrot\qx \ or ek+q < u}-Vrot\qy\^ 
for the two terms, respectively. 

We focus on the self-energy for momenta exactly on 
the Fermi surface, |k| = kp. In this case, the con- 
straints clearly require small q for small uj. For such 
small wavevectors, one may approximate Ck+q ~ vp ■ <i~'r 

/2m* . For a generic point on the Fermi surface (taken 
for simplicity in the quadrant with k^, ky > , for which 
vp makes an angle 9p ^ 0, it/2 to the x axis, the q^ term 
is negligible, and one has Ck+q ~ vpiqx cos 9p+qy sin6'i?). 
Applying the above constraints, one finds that both q^ 
and qy are integrated over a small bounded region in 
which both are 0(w). Thus for such generic points 
on the Fermi surface, we should consider the limit of 
uln)"{c{,Lu — Ck+q) in which all arguments are 0{u!). In 
this limit, inspection of Eqs. (|242I245|I shows that the 
correlator satisfies a scaling form, 

C/(f '(q, LU - ek+q) ~ c.2(A.-2)z^(gj^^ g^/^)^ (247) 

with a well-behaved limit as q ^ 0. Inserting this into 
Eq. H241|) and rescaling the momentum integrals by to, 
one sees that I]"y(ki?, w) ~ uj'^i^^-^') for 9p ^ 0,7r. Since 
Ay > 2 in the stable region of the RFL, this dependence 
is always as weak or weaker than the ordinary uj'^ scatter- 
ing rate due to Coulomb interactions in a Fermi liquid. 

The special cases when 6p ~ 0,7r/2 require separate 
consideration. For these values, the Fermi velocity is 
along one of the principal axes of the square lattice, and 
the linear approximation for Ck+q is inadequate. Taking 
for concreteness 6*^ = 0, we have instead ek+q ~ '^pqx ^" 
g^/2m*, and it is not obviously consistent to neglect the 
qy term. For the first term (involving R"[qy,uj — Ck+q)), 
the constraints reduce to < vpqx ~\- q^/2m* < uj — 
Vrot\qy\- This is approximately solved for small w by < 
qy < {uj-Vrot\qx\)/vp and \qx\ < uj/vp. Thus both q^, qy 
are again bounded and 0{uj), so the above Fermi liquid- 
like scaling applies. 

For the second term (involving R"{qx,uj — ek+q)), the 
constraints reduce to < vpqx + qy/2m* < uj — Vrot\qx\- 



This is solved by taking —uj/vrot < qx < ^/{vrot + vp) 
and Ta&yi{0,-vpqx) < g^/2m* < uj ~ vpq^ - Vrot\qx\- 
Hence for this term, q^ is 0{uj) while qy is 0{^/uj). Since 
with this scaling, \qy \ ^ uj \qx\, one has the significant 
simplification 



Vy{q,uj) 



qx^o 



2w+ v^„tql-uj^ 



for \qy\ > w \qx\. 



Since only Aq appears, in this limit, C/qq"*" '3> ?7, 
i.e. the fluctuations of /Jq arc much stronger than those 
of One may therefore approximate 



(248) 



Oi 



7=0 



4 qxC^ '^rotqx) 



2 

xQ(\uj\ - Vrot\qx\)sg-i^uj5f,Q5^(). (249) 

Inserting this into Eq. (|241|l . one may integrate over qy 
(to yield a constant multiplying y^j^TT and rescale qx — > 
ujqx^ to obtain the result 



Kfi^F.Uj) 



,2A^, 



for( 



TT 

«'2'' 



(250) 



At the end point of the RFL, for which Ay — 2, this 
gives the anomalous lifetime S"(ki?,a;) ~ w'^/^ at these 
special "hot spots" for which the Fermi velocity is parallel 
to one of the principle axes of the square lattice. For 
the conventional Fermi surface believed to apply to the 
cuprates, this corresponds to the points on the Fermi 
surface closest to (7r,0), (0,7r). 

Comparing the scattering rate at these hot spots to 
elsewhere on the Fermi surface, one finds 



Kf;yM 



flUJ 



(251) 



Since the Fermi surface in the cuprates is particularly 
"flat" near (tt, 0), the effective mass would be large, 
TO* ^ TOe, and a significantly enhanced scattering rate 
at such "hot spots" would be expected in the RFL phase. 

In parameter regimes where the RFL phase is unstable 
(via "charge hopping" ) at low temperatures to a conven- 
tional superconductor, the enhanced scattering at the hot 
spots will be significantly suppressed upon cooling below 
the transition temperature. Indeed, the low energy ro- 
ton excitations are gapped out in the superconducting 
phase, and at temperatures well below Tc will not be ap- 
preciably thermally excited. Unimpeded by the rotons, 
the electron lifetime will be greatly enhanced relative to 
that in the normal state, particularly at the hot spots on 
the Fermi surface with tangents along the x or y axes, 
for example at wavevectors near (tt, 0) in the cuprates. 



2. Scattering mediated by "superconducting fluctuations", 
i.e. boson-fermion pair exchange 

A second mechanism of fermion decay is by the "pair- 
ing term" TCa in Eq. (|213|l . We supposed |A| is small. 
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and so consider the first perturbative contribution to the 
quasiparticle lifetime, at O(A^). In general, for d-wave 
pairing, this takes the form 



again with possible analytic and sub-dominant correc- 
tions. Upon analytic continuation to obtain the lifetime, 
only non-analytic terms contribute, and we expect 



2 

E ^^^^ 



sign(cj). 



(257) 



qf2„ 



-i{^n - ^n) + Cq-k' 



where 



(r,r)e 



zq-r — lOjj T 



with G'^j'(r,r) from Eq. ^U^. 

Our uncertainties in the details of the Cooper pair 
propagator (originating from ambiguities in the string 
geometry) do not allow a thorough calculation of the re- 
sulting self energy. However, as we will now show, we can 
obtain a rough understanding of its scaling properties by 
some simple approximations, which we believe do not sig- 
nificantly effect the results. In particular, we will assume 
that the lifetime is controlled by the small wavevector 
\q\ <C TT portion of the Cooper pair propagator. In this 
regime, we expect G'^J(q, ri„) « ^^^(q, r2„), the Fourier 
transform of the local Cooper pair propagator studied in 
depth in Sec. IIII CI Making this approximation, one has 
l^l for k near the quasiparticle Fermi surface, and 

one may write 



for k on the Fermi surface. This result can be verified by 
more detailed calculations using the spectral representa- 
on of G'^^M' Moreover, we expect that for kgT ^ uj, oj 
can be replaced by ksT in this formula. 

As expected, for sufficiently large Ac, this lifetime 
vanishes rapidly at low energies, and in particular for 
(253) Ac > 3/2, this contribution is sub-dominant to the usual 
Fermi liquid one. However, in the regime with Ac < A„, 
this is never the case (Ac < l/^/2), at least within our 
simple model without dramatic corrections to the RL ex- 
ponents. Indeed, for Ac < 1, as supposed herein, this 
"scattering rate" is much larger than the quasiparticle 
energy at low frequency. Taken literally, this implies in- 
creasingly incoherent behavior away from the nodes as 
temperature is lowered. Near the nodes, the amplitude 
of this strong scattering contribution vanishes rapidly, 
similar to the idea of "cold spots" proposed by loffe and 
MiUisia 



I]s/(k,cj„) 



-i(ri„ 



-(254) 



with Ak = |A|(cosfca; — cos ky). Note that this form im- 
mediately implies that scattering due to this mechanism 
is strongly suppressed upon approaching the nodal re- 
gions of the Fermi surface. 

A detailed analysis is now possible based on a spec- 
tral representation of G"^p, as in Sec. IVII B itf Unlike 
the above case, however (since the power law decay of 
G"^P(r, r) is approximately isotropic) a much simpler scal- 
ing analysis sufhces to obtain the qualitative behavior of 
the lifetime. In particular, the scaling form of Eq. H127|) 
implies that 



G^P(q,r!„)~ la. 



^(q/f^n), 



(255) 



up to non-singular additive corrections, for small \q\ and 
ri„. Furthermore, for small q and k on the Fermi surface, 
one may write Cq^k ~ Vpqw +q'j^/2m, where q\\ , q± are the 
components of q parallel and perpendicular to the local 
Fermi velocity at k, and Vp and m are the magnitude of 
the local Fermi velocity and effective mass. The singu- 
larity in the integrand in Eq. H254|) is then cut off by the 
external frequency w„, and rescaling q qr2„, one sees 
that the effective mass term is negligible, and obtains by 
power counting. 



i;s/(k, w„) ~ |Ak| 



|2Ac-l 



VIII. DISCUSSION 



We close with a discussion of some theoretical issues 
concerning the vortex-quasiparticle formulation of inter- 
acting electrons that we have been employing through- 
g^ j^ - ' out. In particular, we contrast our approach with the 
earlier Z2 formulation^ and mention the connection with 
more standard and more microscopic formulations of cor- 
related electrons. We then address the possible relevance 
of the Roton Fermi Liquid phase to the cuprate phase di- 
agram and the associated experimental phenomenology. 



A. Theoretical Issues 

Since the discovery of the cuprate superconductors, 
the many attempts to reformulate theories of two- 
dimensional strongly correlated electron in terms of new 
collective or composite degrees of freedom, have been fu- 
eled on the one hand by related theoretical successes and 
on the other by cuprate phenomenology. The remark- 
able and successful development of bosonizationi^ both 
as a reformulation of interacting one-dimensional elec- 
tron models in terms of bosonic fields and as a means to 
extract qualitatively new physics outside of the Fermi 
liquid paradigm, played an important role in a num- 
ber of early theoretical approaches to the cupratesi^. 
The equally impressive successes of the composite bo- 
son and composite fermion approaches in the fractional 
quantum Hall effectSS were also influential in early high 
Tc theories, most notably the anyon theories2i. Gauge 
(256) theories of the Heisenberg and t — J models were no- 
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table early attempts to reformulate 2d interacting elec- 
trons in terms of "spin-charge" separated variables^SiS^ 
- electrically neutral fermionic spin one-half "spinous" 
and charge e bosonic holons - and were motivated both 
by analogies with Id bosonization and by resonating 
valence bond ideasi^. More recent approaches to 2d 
spin-charge separation have highlighted the connections 
with superconductivity"*"^, by developing a formulation 
in terms of vortices, Bogoliubov quasiparticles and plas- 
mons: the three basic collective excitations of a 2d su- 
perconductor. These latter theories were primarily at- 
tempting to access the pseudo-gap regime by approaching 
from the superconducting phased - focusing on low en- 
ergy physics where appreciable pairing correlations were 
manifest. 



1. Vortex-fermion formulation 

In this paper, although we are employing a formulation 
with the same field content - hc/2e vortices, fermionic 
quasiparticles and collective plasmons - we are advocat- 
ing a rather different philosophy. In particular, we wish 
to use these same fields to describe higher energy physics 
in regimes where the physics is decidedly non-BCS-like 
even at short distances, most importantly the cuprate 
normal state near optimal doping. The philosophy of this 
approach is very similar to that of the Z2 gauge theory 
proposal of a fractionalized under-doped normal state. 
Indeed, as demonstrated in Sec. Ill Cl and Adds. IbIcI our 
C/(l) formulation is completely (unitarily) equivalent to a 
Z2 gauge theory. However, because the unitary transfor- 
mation relating the two formulations is non-trivial, the 
U(l) formulation is (much) more convenient for the types 
of manipulations and approximations we employ here, 
largely because the C/(l) gauge fields are continuous vari- 
ables. The price paid for the use of the U{1) formulation 
is that the "spinon pairing term" of the Z2 gauge theory 
appears non-local in the U{1) Hamiltonian. 

Fortunately, this non-locality and its consequences 
are readily understood. In particular, although the 
C/(l) vortex-quasiparticle Hamiltonian is microscopically 
equivalent to a Z2 gauge theory, it is also equivalent (as 
described in Add. ID|| to a theory of electronic (i.e. charge 
e) quasiparticles coupled to charge 2e bosons (with the 
latter described in dual vortex variables). The "two- 
fluid" point of view of this vortex-electron formulation is 
convenient for understanding the qualitative properties of 
the RFL phase and its descendents, although it should be 
emphasized that the current-current couplings (embod- 
ied by the gauge fields of the U{1) formulation) between 
the two fluids are not expected to be weak. In the vortex- 
electron language, the non-local term simply represents 
coherent "Josephson coupling" of electron pairs and the 
bosons. The non-locality of the pairing term arises sim- 
ply from the non-local representation of boson operators 
in the usual U{1) boson-vortex duality. The crucial fea- 
ture which allows us to handle the pairing term despite 



its non-locality is the strength of vortex fluctuations in 
the bosonic Roton Liquid, which persists even with these 
strong current-current couplings. The resulting power- 
law decay of bosonic pair-field correlations (ODQLRO) 
opens up a regime, corresponding to the RFL phase, in 
which the pairing term is irrelevant and can be treated 
perturbatively. We note in passing that, although it is 
not relevant for the cuprates, the above discussion makes 
clear that a RL phase (for which fermions need never be 
introduced) should be possible in a purely bosonic model, 
which would be interesting in and of itself. 

A second important feature of the present formulation 
is the retention of the lattice scale structure. Appropriate 
to the cuprates, we have carefully defined the theory on 
a 2d square lattice which we wish to identify with the mi- 
croscopic copper lattice. Since our theory obviously does 
not similarly retain physics on atomic energy scales, our 
starting "bare" Hamiltonian should be viewed as a low 
energy but not spatially coarse grained effective theory 
(or at most an effective theory coarse-grained spatially 
only insofar as to remove e.g. the O p-orbitals). It is im- 
portant to emphasize that the very existence of the Ro- 
ton Fermi liquid phase requires the presence of a square 
lattice. In the RFL ground state there is an infinite set 
of dynamically generated symmetries, corresponding to 
a conserved number of vortices on every row and column 
of the lattice. If we study the same model on a different 
lattice, say the triangular lattice, the quantum ground 
state analogous to the RFL phase (obtained, again, by 
taking the roton hopping amplitude Kr 00) is highly 
unstable, being destroyed by the presence of an arbitrar- 
ily small single vortex hopping amplitude. 

More generally, the importance of employing the dual 
vortex-quasiparticle field theory reformulation cannot be 
overemphasized. The novel and unusual RFL phase 
emerges quite simply when one takes a large amplitude 
for the roton hopping term, and the fixed point theory is 
also quite simple consisting of a Fermi sea of quasiparti- 
cles minimally coupled to an "electric" field with Gaus- 
sian dynamics. It is very difficult to see how one could 
adequately describe such a phase working with bare elec- 
tron operators. Indeed, we do not yet understand the 
simpler task of constructing a microscopic (undualized) 
boson model which enters the RL phase, even without the 
complications of fermions. However, previous experience 
with dual vortex formulations for bosonic and electronic 
systems strongly argues that the RL and RFL theories 
properly describe physically accessible (albeit microscop- 
ically unknown) models. 

An unfortunate drawback of the present formulation, 
however, is the apparent lack of direct connection be- 
tween the starting lattice Hamiltonian and any micro- 
scopic electron model. In particular, it is presently un- 
clear what microscopic electron physics could be respon- 
sible for generating such a large roton hopping term. In 
the absence of a microscopic foundation, one must re- 
sort to developing phenomenological implications of the 
Roton Fermi Liquid phase and comparing with the ex- 
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perimentally observed cuprate phenomenology. We take 
a preliminary look in this direction in Sec. IVlll Bl below. 



2. Related approaches 

It is instructive to contrast the RFL phase with ear- 
lier notions of spin-charge separation. Anderson's origi- 
nal picture of the cuprate normal state at optimal dop- 
ing consisted of gas of spinous with a Fermi surface co- 
existing and somehow weakly interacting with a gas of 
holons^'^'^ In the early gauge theory implementations of 
this picture, both the spinous and holons carried a U{1) 
(or SU{2)) gauge charg o^^i^^ . Mean-field phase diagrams 
were obtained by pairing the spinous and/or condens- 
ing the holons, and a pseudo-gap regime with d-wave 
paired spinous yet uncondensed holons was predicted^^iSi 
- several years before experiment (of course the d-wave 
nature of the superconducting state was also predicted 
by other approaches^). Within this approach, the nor- 
mal state at optimal doping was viewed as an incoher- 
ent gas of uncondensed holons strongly interacting with 
unpaired spinons. Some efforts were made to use the 
neglected gauge field fluctuations to stop the holons con- 
densing at inappropriately high energy scales, with some 
succes9iS4. A serious worry is that these gauge fluctu- 
ations, ignored at the mean-field level, will almost cer- 
tainly drive confinement?^ (gluing the spinons and the 
holons together) and thus at low temperatures invali- 
date the assumed stability of the initial mean-field saddle 
point. 

A few early papers emphasized that spinon pairing 
would break the continuous U{1) (or SU{2)) gauge 
symmetry down to a discrete Z2 subgroupi^, and 
this might be a way to avoid confinement. These 
ideas were later considerably amplified^, and the sta- 
bility of genuinely spin-charge separated ground states 
established^'^"^^. Such ground states are exotic electrical 
insulators which support fractionalized excitations carry- 
ing separately the spin and charge of the electron. Each 
fragment carries a Z2 gauge charge, but in contrast to 
the U{1) gauge theory saddle points, the Z2 spinons are 
electrically neutral and do not contribute additively to 
the resistance. Such fractionalized insulators necessar- 
ily support an additional Z2 vortex-like excitation: the 
vison^. 

Within the present formulation, these Z2 fractionalized 
insulators can be readily accessed by condensing pairs of 
vortices'^, rather than condensing rotons. Since vortex- 
pairs only have statistical interactions with charged par- 
ticles and not the spinons, the resulting pair-vortex con- 
densate is an electrical insulator with deconfined spinon, 
"chargon" and vison excitationsiii^ - dramatically differ- 
ent from the Roton Fermi Liquid. Recent experiments on 
very under-doped cuprate samples have failed to find evi- 
dence of a gapped visor^l^^ , suggesting that the pseudo- 
gap phase is not fractionalized. But a negative result in 
these vison detection experiments does not preclude the 



RFL phase, which supports mobile single vortices even at 
very low temperatures. A different approach supposing 
unbound and mobile single vortices is the QED3 theory 
of Ref. Ti. 

At the most basic level, a roton is a small pattern of 
swirling electrical currents. Much recent attention has 
focused on the possibility of a phase in the under-doped 
cuprates with non- vanishing orbital currents'^*', counter 
rotating about elementary plaquettes on the two sub- 
lattices of the square lattice. Such a phase was initially 
encountered as a mean-field state within a slave-Fermion 
gauge theory approach^! - the so-called "staggered-flux 
phase" , but has been resurrected as the "d-density wave" 
ordered state of a Fermi liquidi2&. In either case, such a 
phase within the present formulation would be described 
as a "vortcx-antivortex lattice" - a checkerboard configu- 
ration on the plaquettes of the 2d square lattice. Ground 
states with short-ranged "orbital antiferromagnetic" or- 
der have also been suggested recently??. Surprisingly, 
the Roton liquid phase also has appreciable short-ranged 
orbital order. A "snapshot view" of the orbital current 
correlations in the RL phase can be obtained by examin- 
ing the vortex-density structure function. For simplicity, 
consider the charge sector of the RFL theory (the RL) in 
the limit of large plasmon velocity, uq — > 00, which sup- 
presses the longitudinal charge density fluctuations. The 
transverse electrical currents are then described by the 
Hamiltonian Hrot in Eq- (|89|1 with aj 0. The vortex- 
density structure function, Snn, can be readily obtained 
from this Gaussian theory: 

5^^.(k) ^ (|iV(k)n = iy^|/C.(k)/C,(k)|v/^, 

(258) 

with \ICj{k)\ = 2|sin(fcj/2)| and IC^{k) = l^jWP 
as before. The vortex structure function is analytic ex- 
cept on the kx = and ky = axes and is maxi- 
mum at k = (tt, tt), indicative of short-ranged orbital- 
antiferromagnetism. As such, the RFL/RL phase can 
perhaps be viewed as a quantum-melted staggered-flux 
(or vortex-antivortex) state, with only residual short- 
ranged orbital current correlations - the rapid motion 
of the rotons being responsible for the melting. If the 
fcrmions are paired, the amplitude for the basic roton 
hopping process which generates the structure function 
above vanishes upon approaching half filling, as is ap- 
parent from Eq. (|84|l . A preliminary analysis^'^ suggests 
that the further neighbor roton hopping processes which 
do survive at half-filling, lead to a vortex structure func- 
tion which vanishes as |k — (7r,7r)| for wavevectors near 
(7r,7r). 



B. Cuprates and the RFL phase? 

In applying BCS theory to low temperature supercon- 
ductors, one implicitly assumes that the normal state 
above Tc is adequately described by Fermi liquid theory. 
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Within a modern renormalization group viewpoint'^ this 
is tantamount to presuming that the effective Hamilto- 
nian vaHd below atomic energy scales (of say lOeV) is 
sufficiently "close" (in an abstract space of Hamiltonians) 
to the Fermi liquid "fixed point" Hamiltonian (actually 
an invariant or "fixed" manifold of Hamiltonians char- 
acterized by the marginal Fermi liquid parameters). In 
practical terms, "close" means that under a renormaliza- 
tion group transformation which scales down in energies, 
the renormalized Hamiltonian arrives at the Fermi liquid 
fixed point on energy scales which are still well above T^. 
BCS theory then describes the universal crossover flow 
between the Fermi liquid fixed point (which is marginally 
unstable to an attractive interaction in the Cooper chan- 
nel) and the superconducting fixed point which charac- 
terizes the universal low energy properties of the super- 
conducting state well below Tc- Four orders of magni- 
tude between f OeV and Tc gives the RG flows plenty of 
"time" to accomplish the first step to the Fermi liquid 
fixed point, and is the ultimate reason behind the amaz- 
ing quantitative success of BCS theory. 



1. Assumptions underlying the RFL approach 

In what follows, our working hypothesis is that the 
effective electron Hamiltonian on the lOeV scale appro- 
priate for the 2d copper-oxygen planes at doping lev- 
els within and nearby the superconducting "dome" , is 
"close" to the Roton Fermi liquid "fixed manifold" of 
Hamiltonians (characterized by Fermi/Bose liquid pa- 
rameters for the quasiparticles/rotons). More specifi- 
cally, we presume that when renormalized down to the 
scale of say one-half of an eV, the effective Hamiltonian 
can be well approximated by a Hamiltonian on the RFL 
fixed manifold - up to small perturbations. The impor- 
tant small perturbations are those that are relevant over 
appreciable portions of the fixed manifold. As established 
in Sections IV and V, there are three such perturbations: 
(i) Single vortex hopping, (ii) "charge" hopping and (iii) 
attractive quasiparticle interactions in the Cooper chan- 
nel. When relevant, these three processes destabilize the 
RFL phase and cause the RG flows to cross over to dif- 
ferent fixed points which determine the asymptotic low 
temperature behavior. From our calculations, we find 
that at least one, and often more of these three pertur- 
bations is always relevant, regardless of the roton liquid 
parameters. Hence the RFL should be regarded as a criti- 
cal, and usually multi-critical phase, rather than a stable 
one. For these three processes, the resulting quantum 
ground states are, respectively: (i) A ("confined" and) 
conventional Fermi liquid phase, (ii) a conventional su- 
perconducting phase with singlet d^2_y2 pairing and gap- 
less nodal Bogoliubov quasiparticles and (iii) a "rotonic 
superconductor" with the properties of a conventional 
superconductor but coexisting gapless roton excitations. 
The rotonic superconductor, as described in Sec. I VI B2l 
has further potential instabilities driven by either single 



vortex hopping and "charge" hopping/explicit pairing. 
Both perturbations, if relevant, will generate an energy 
gap for the rotons. It seems likely that the domains of rel- 
evance of these two perturbations overlap, so that there 
no regime of true stability of the rotonic superconductor. 
The true ground state of the system in this regime is then 
a conventional superconductor, and its "rotonic" nature 
is evidenced only as an intermediate energy crossover. 

2. Effective parameters 

In order to construct a phase diagram within this sce- 
nario, it is necessary to specify the various parameters 
(eg. Bose and Fermi liquid parameters) of the RFL model 
as a function of the doping level, x. Because the RFL 
is multi-critical, we cannot rely upon "universality" to 
validate ad-hoc requirements of smallness of perturba- 
tions, as might be the case e.g. for renormalized per- 
turbations around a stable fixed point. Instead, we will 
make some assumptions based (partly) on physics. First, 
our main assumption is the basic validity at high energies 
of the roton dynamics, and of Fermi liquid-like quasipar- 
ticles. Second, we assume that superconductivity is never 
strong, i.e. always occurs below the rotonic/Fermi liquid 
energy scales. Mathematically, these two assumptions 
are encompassed in the inequalities 

Kr,Uv,ts tv :S> tc,Aj. (259) 

Reading from left to right, this corresponds to the first 
and second assumptions above. In practice, we can at 
best hope for a factor of a few between e.g. ty and Hr, so 
the ":^" symbols above should not be taken too strongly. 
Although it is not important to our discussion, it is natu- 
ral to assume that vq Kr.,Uy, tg (since Coulombic ener- 
gies at the lattice scale are comparable to the electronic 
bandwidth). A third assumption, which is not needed 
for consistency of the approach, but seems desirable em- 
pirically, is that the fermion dynamics is primarily by 
"spinon" rather than electron hopping, tg 3> ie- Many 
of the parameters of the RFL phase can be fixed empir- 
ically from the observed behavior of the cuprates on or 
above the eV scale. For example, the fc— space location 
of the quasiparticle Fermi surface can be chosen to co- 
incide with the electron Fermi surface as measured via 
ARPES^. The value of other parameters, such as the 
bare velocity vq which appears in Ca and sets the scale of 
the plasmon velocity, can be roughly estimated from the 
basic electronic energy scales, and in any case does not 
greatly effect the relevance/irrelevance of the three im- 
portant perturbations. For our basic lattice Hamiltonian 
introduced in Sec. ^ the two most important param- 
eters characterizing the RFL phase are the vortex core 
energy, and the roton hopping strength, k^. Indeed, 
as shown in Sec. IIVI the scaling dimensions which deter- 
mine the relevance of the vortex and "charge" hopping 
perturbations at the RFL fixed point, denoted A^, and 
Ac respectively, depended sensitively on the ratio u^/Kr- 
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FIG. 2: Proposed values for the vortex and "charge" hopping 
scahng dimensions, A„(j:) and Ac{x), as a function of the 
doping X. These hopping perturbations are relevant at the 
RFL fixed point when their A < 2. Within the doping range, 
xi < X < X2, the "charge" hopping is more relevant, and the 
RFL phase is unstable to superconductivity. 



For example, ignoring renormalizations from the gapless 
fermions we found. 



1 

2AI 



(260) 



with v'^ — + jUyKr ^ Vq + -^UjjKr- 

Rather than specifying the doping dependence of 
Uy, Kr, vq and the Bose/Fermi liquid parameters, though, 
it is simpler and suffices for our purposes to specify the 
final X— dependence of A^, and Ac. Shown in Figure 2 is 
a proposed form for Ay/^i^), primarily chosen to fit the 
gross features of the cuprate phase diagram. For exam- 
ple, since the RFL phase is strongly unstable to the Fermi 
liquid phase when A„ <C 2, we have taken A^ decreasing 
to small values in the strongly over-doped limit. On the 
other hand, to account for the observed non-Fermi liq- 
uid behavior in the normal state near optimal dopingi^, 
requires that we take A^ > 2 in that regime. On the 
under-doped side of the dome we can use the observed 
linear x-dependence of the superfluid density to guess the 
behavior of A^ for small a;. In this regime the (renormal- 
ized) vortex core energy in the superconducting state is 
presumably tracking the transition temperature, varying 
linearly with x. It seems plausible that the "bare" vor- 
tex core energy Uy in the lattice Hamiltonian, while per- 
haps significantly larger, also tracks this x— dependence. 
This implies A^, ^ ,Ju^ ~ ^/x as depicted in the Figure. 
Moreover, to recover (conventional) insulating behavior 
when a; — > 0, requires that vortex condensation (rather 
than charge hopping) be more relevant in this limit, i.e. 
A^ < Ac (see below). 



3. Phase diagram 

Under the above assumptions, we now discuss in some 
detail the resulting phase diagram and predicted behav- 
iors. Consider first the ground states upon varying x. In 



the extreme over-doped limit with A„ <C 2, vortex hop- 
ping will be a strongly relevant perturbation at the RFL 
fixed point. The vortices will condense at T = 0, lead- 
ing to a conventional Fermi liquid ground state. Upon 
decreasing x there comes a special doping value {x2 in 
Fig. 2) where Ac becomes smaller than A^,. Ki x = X2 
the RFL phase is unstable to both vortex and "charge" 
hopping processes, since both A„ = Ac = l/\/2 < 2. 
It seems reasonable to assume that in this situation, the 
more strongly relevant process ultimately dominates at 
low energies. This implies that X2 demarcates the bound- 
ary between a Fermi liquid and a superconducting ground 
state, as illustrated in Figure 3. Upon further decreasing 
a;, the "charge" hopping becomes even more strongly rel- 
evant, tending to increase Tc until it reaches a maximum 
at "optimal doping" , denoted Xopt in Figures 2 and 3. As 
one further decreases x, Tc should start decreasing. But 
as shown in Sec. VA, with decreasing vortex core energy, 
Ml, ~ X, the enhanced vortex density fluctuations gener- 
ate an increasing antiferromagnetic exchange interaction 
J ~ t'^/uy. This attractive interaction between fermions 
will mediate quasiparticle pairing, with a pairing energy 
scale growing rapidly as a; — > 0. It is natural to asso- 
ciate this "quasiparticle pairing" temperature scale with 
the crossover temperature into the pseudo-gap regime—, 
shown as T* in Figure 3. As discussed in Sec. IVI B 21 
under the assumption that quasiparticle kinetic energy 
arises primarily through spinon hopping, tg ^ t^, the 
superfluid density associated with the quasiparticle pair- 
ing is small, so that potential true superconductivity as 
a consequence of this pairing is suppressed to a low or 
zero temperature. 

Since A„ < Ac for doping levels with a; < a;i, vor- 
tex hopping should again dominate over "charge" hop- 
ping. This is the same condition which we argued leads 
to the Fermi liquid state for a; > a;2 above. However, 
the physics for small x is more complex, owing to the 
strong antiferromagnetic interactions and proximity to 
the commensurate filling a; = at which antiferromag- 
netic order is probable. In the U{1) vortex-quasiparticle 
formalism of this paper, this difference arises from the 
freedom to choose (as Aj 0) the fermion density to 
minimize the total (free) energy of the system. In the 
majority of this paper we have taken nj — po, in order to 
minimize the vortex kinetic energy. However, if antiferro- 
magnetic quasiparticle interactions are large, and a; <C 1, 
another possibility arises. To optimally benefit from the 
antiferromagnetic interactions, one may instead choose 
the fermion density commensurate, Uf = 1, for which 
the Fermi surface is optimally nested and the quasiparti- 
cles can become fully gapped, gaining the maximal "con- 
densation energy" from antiferromagnetic ordering. The 
cost of this choice is some loss of vortex kinetic energy as 
the vortex motion becomes somewhat frustrated by the 
resulting dual "ffux" Tra:. Although we assume vortex en- 
ergy scales are large, this fiux itself is small for small x, 
so eventually as a; — s- this rise in kinetic energy becomes 
smaller than the lowering of fermionic energy due to an- 
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FIG. 3: Schematic phase diagram that follows from the dop- 
ing dependence of A„,c shown in Fig. 2. With increasing 
doping X, the ground state evolves from a charge ordered in- 
sulator (I), into a superconductor (dSC) and then back into 
a Fermi liquid (FL). The normal state behavior near optimal 
doping is controlled by the Roton Fermi liquid (RFL) ground 
state. Below T*, the quasiparticles pair. A quantum phase 
transition between the RFL and FL ground states occurs at 
Xc - but is preempted by superconductivity. At x = Xc, a. 
normal state resistance linear in temperature is predicted. 



tiferromagnetism and such a choice becomes favorable. 
The ultimate physics of the remaining vortices in this 
limit is difficult to analyze reliably, but it seems evident 
that for < Ac, the more relevant vortex hopping will 
lead to an insulating state. Since the quasiparticles are 
then gapped, the remaining uncompensated flux ttx is ex- 
pected to lead to incommensurate charge ordering when 
the vortices condense - the dual analog of the Abrikosov 
lattice. 



4. The RFL Normal state 

Having established the phase diagram which follows 
from the assumed doing dependence of A^, (x) , we turn to 
a discussion of the normal state properties above Tc. Un- 
der our working assumption, in the energy range above 
Tc we can ignore the small (but ultimately relevant) 
"charge" hopping perturbation, and use the RFL fixed 
point Hamiltonian to describe the physics of the normal 
state. First consider the special doping value, x = Xc, 
where A^, = 2 on the over-doped side of the supercon- 
ducting dome (see Figures 2 and 3). For x < Xc the 
vortex hopping strength decreases upon scaling down 
in energy, so that it can be treated perturbatively. As 
shown in Sec. VIA, at second order in the vortex hop- 
ping, the electrical conductivity is additive in the roton 
and quasiparticle conductivities. With impurities present 
the quasiparticle conductivity, ^ will saturate to low 
temperatures, whereas Urot ^ T^'^ (with 7 = 2A^ — 3) 



diverges as T ^ 0. The low temperature normal state 
electrical resistance is thus predicted to behave as a power 
law. 



R{T) - tlT'^, 



(261) 



with a vanishing residual resistivity. Moreover, right at 
X = Xc, since 7 = 1 a linear temperature dependence is 
predicted. Notably, this transport behavior is a due to 
the presence of the quasi-condensate in the RFL phase, 
and as such is completely independent of the single par- 
ticle scattering rate. 

The Hall conductivity, however, will be largely deter- 
mined by the fermionic quasiparticle contribution, as de- 
tailed in Sec. IVIII Specifically, the cotangent of the 
Hall angle, cot(0/f) = Pxx/Pxy, was found to vary as 
cot(eH) ~ l/iT'^T]) in the RFL phase, with tJ^ the 
spinon momentum relaxation rate. Moreover, at the "hot 
spots" on the Fermi surface, tJ^ ~ T'^'+s^ due to scat- 
tering off the gapless rotons. Under the assumption that 
these same processes dominate the temperature depen- 
dence of the quasiparticle transport scattering rate, we 
deduce that. 



cot(e^f) 



1 



(262) 



and right at quadratic dependence cot(O^f) ~ 

T^. In striking contrast to conventional Drude theory 
which predicts cot(8H) ^ ^ R (with r the elec- 
tron's momentum relaxation time), in the RFL phase 
the cotangent of the Hall angle varies with a different 
power of temperature than for the electrical resistance, 
R ^ T'' . This non-Drude behavior is consistent with the 
electrical transport generally observed in the optimally 
doped cuprates2S*^Si2i, where i? ~ T and cot(eff) ~ T^. 

Consider next the thermal conductivity k near optimal 
doping within the RFL normal state. One of the most im- 
portant defining characteristics of a conventional Fermi 
liquid is the Wiedemann-Franz law - the universal low 
temperature ratio of thermal and electrical conductivi- 
ties, L = n/aT. In a Fermi liquid, electron-like Landau 
quasiparticles carry both the conserved charge and the 
heat, and since the energy of the individual quasiparti- 
cles becomes conserved as T 0, the Lorenz ratio is 
universal, L^^ = Lq — Tr^kg/3e^. In contrast, the elec- 
trical conductivity is infinite in a superconductor, but 
the condensate is ineffective at carrying heat so that the 
Lorenz ratio vanishes, L*"^ = 0. Within the RFL phase, 
heat can also be transported by the single fermion exci- 
tations, with a contribution to the thermal conductivity 
linear in temperature: Kg = L^a^T . At low enough tem- 
peratures this will dominate over the phonon contribu- 
tion, Kphon ~ T^- But the roton excitations, which have 
a quasi-one dimensional dispersion at low energies, will 
presumably also contribute a linear temperature depen- 
dence, Krot ^ T. Thus, the total thermal conductivity in 
the RFL phase is expected to vanish linearly in temper- 
ature, K T . But since the roton contribution to the 
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electrical conductivity diverges as T ^ 0, the RFL phase 
is predicted to have a vanishing Lorenz ratio: 

^BFL ^ Jl. _ (263) 

The quasi-condensate in the RFL phase is much more 
effective at transporting charge than heat, much as in a 
superconductor. Electron doped cuprates near optimal 
doping, when placed in a strong magnetic field to quench 
the superconductivity, do exhibit a small Lorenz ratio 
at low temperatures^^, L « Lo/5. But extracting the 
zero field Lorenz ratio is problematic, since above Tc the 
phonon contribution to k is non-negligible. 

It is instructive to consider the electrical resistance 
also in the under-doped regime, particularly upon cool- 
ing below the fermion pairing temperature. Above this 
crossover line the predicted electrical resistance varies 
with a power of temperature, R ~ T'^ . Since, as we 
assume, the fermion's kinetic energy comes primarily in 
the form of spinon hopping, t k tg ^ te, the resulting 
superfluid density is however very small (see Sec. lVIB^ . 
and phase coherent superconductivity does not result, at 
least not in this temperature range. Nevertheless, one 
would expect a dramatic increase in the fermion con- 
ductivity, , upon cooling through T* - much as seen 
in superconducting thin films upon cooling through the 
materials bulk transition temperature. Since the conduc- 
tivity is additive in the roton and spinon contributions, 
R{T)^^ ^ arot{T) + <y^ {T), a large and rapid increase 
in <J-f{T) should be detectable as a drop in the electrical 
resistance relative to the "critical" power law form, i.e. 

R{T) 1 

^ l + cT^afiT)' ^^^"^^ 

This behavior is generally consistent with that seen in 
the under-doped cuprates^i^, provided we take 7 « 1. 

5. Entering the superconductor 

We finally discuss the predicted behavior upon cool- 
ing from the RFL normal state into the superconducting 
phase. The main change occurs in the spectrum of ro- 
ton excitations, which become gapped inside the super- 
conductor. The roton gap, A^ot, should be manifest in 
optical measurements, since the optical conductivity will 
drop rapidly for low frequencies, uj < Aj-ot- As in BCS 
theory, the ratio of the (roton) gap to the superconduct- 
ing transition temperature 2Arot/Tc should be of order 
one. This ratio is determined by the RG crossover flow 
between the RFL and superconducting fixed points. It 
seems likely that this ratio will be non-universal, depend- 
ing on the marginal Bose liquid parameters characteriz- 
ing the RFL fixed point (this is distinct from changes 
in this ratio in strong coupling Eliashberg theories, since 
here variations in the ratio are due to marginal parame- 
ters of the RFL fixed manifold even at arbitrarily weak 



coupling). Nevertheless, it would be instructive to com- 
pute this ratio for the simple near-neighbor RFL model 
we have been studying throughout, and to analyze the 
behavior of the optical conductivity above the gap. 

Another important consequence of gapped rotons be- 
low Tc, is that the electron lifetime should rapidly in- 
crease upon cooling into the superconducting phase. 
With reduced scattering from the rotons, the ARPES 
line width should narrow, most dramatically near the 
normal state "hot spots" (with tangents parallel to the x 
or y axes). This behavior is consistent with the ARPES 
data in the cuprates^, which upon cooling into the su- 
perconductor does show a significant narrowing of the 
quasiparticle peak, particularly so at the Fermi surface 
crossing near momentum (7r,0). 

Despite these preliminary encouraging similarities be- 
tween various properties of the RFL phase and the 
cuprate phenomenology, much more work is certainly 
needed before one can establish whether this exotic non- 
Fermi liquid ground state might actually underlie the 
physics of the high temperature superconductors. We 
have already emphasized the strengths of this proposi- 
tion, but there are, of course, some experimental features 
which seem challenging to explain from this point of view. 
The "quasiparticle charge", i.e. temperature derivative of 
the superfluid density dKs/dT\rp^Q appears, based upon 
a small number of experimental data points, to be large 
and roughly independent of doping x, in apparent conflict 
with the RFL prediction. The linear temperature depen- 
dence of the electron lifetime 1/t/ ~ /csT/Ti observed 
for nodal quasiparticles near optimal doping in ARPES 
also does not seem natural in the RFL. However, it seems 
likely that it may be possible to explain a small number 
of such deviations from the most naive RFL predictions 
by more detailed considerations. Further investigations 
of the RFL proposal should confront other experimen- 
tal probes, such as interlayer transport and tunneling. 
Towards this end, it will be necessary to generalize the 
present approach to three-dimensions. More detailed pre- 
dictions, such as for the optical conductivity upon enter- 
ing the superconductor and the thermal Hall effect in the 
RFL normal state, might also be helpful in this regard. 
It would of course be most appealing to identify a new 
"smoking gun" experiment for the Roton Fermi liquid 
state, analogous to the vison-trapping experiment^ for 
detecting 2d spin-charge separation. However, given the 
critical nature of the RFL, with copious gapless excita- 
tions with varieties of quantum numbers, finding such an 
incontrovertible experimental signature may be difficult. 
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APPENDIX A: FERMI LIQUID PHASE 

We expect the Fermi Liquid phase to occur upon com- 
plete proliferation and unbinding of vortices. To ob- 
tain the Fermi liquid in our formulation, we therefore in 
this appendix consider the limit of large vortex hopping 
t-u ^ 00 and small vortex energy uq — > 0. We have pre- 
viously demonstrated the equivalence of the U(l) gauge 
theory to a Z2 gauge theory, and it is this latter formu- 
lation which is most convenient in this limit. 

To observe the Fermi liquid, we analyze the Z2 gauge 
theory in its Hamiltonian form. Although this limit 
is very straightforward, based on previous work on Z2 
gauge theory, it is instructive to go through it in some de- 
tail here, in order to observe the effects of the spinon pair- 
ing term. The Hamiltonian density can be separated into 
pure gauge, charge, and spin parts, H = Hh+'H^'^ +Wf^, 
with 

Hft = -/i^a|(r), (Al) 

i 

where aj is the usual Pauli matrix in the space of states 
on link in the j direction coming from site r (and hence 
anticommutes with a-|). In the charge sector, 

Wf^ = -2tc ^ 0-| (r) C0s((/)r+xj - <^r) + Ucllrlrir - 1], 

(A2) 

where rir is the number operator conjugate to (f)r, satis- 
fying [0r,'T^r'] = iSrr' ■ For simplicity, in the spin sector 
we consider a local s-wave pair field instead of a d-wave 
one. This is not essential, but simplifies the presentation 
and still addresses the essential issue of the relevance of 
spinon pairing in the Fermi liquid. Hence, 

Wf= = -t.a|(r)[/J+-^.,/r, + h.c] + eoflh, 

+A(/rT/ri + h.c.) + UfiflUj. (A3) 

Note that we have added a local on-site energy eo and 
interaction Uf, allowed in general by symmetry. The 
Hamiltonian commutes with the gauge generators, 

Gr = (-1)"-+^-^- [] (r)^J(r - Xj). (A4) 

j 

We require Gr = 1 to enforce gauge invariance. 

In the analysis of the large vortex hopping limit above, 

we obtained the action for a Z2 gauge theory with zero 
kinetic term. This corresponds in the Hamiltonian to 



large h (in particular we will take h S> tc.i.s- In this 
limit crj ~ 1, and the gauge constraint becomes 

i.e. requiring an even number of bosons and fcrmions on 
each site. Further, for large h, the chargon and spinon 
hopping terms are strongly suppressed, and can be con- 
sidered perturbatively. At zeroth order in tc,ts, then, 
the charge and spin sectors are decoupled at each site 
and decoupled also from one another except by the gauge 
constraint. In the charge sector, for Uc > 0, it is energet- 
ically favorable to have only either zero or one chargon 
per site rir = 0, 1. If rir = 0, then we must have either 
zero or two spinous per site. Note than in this subspace, 
even for small A, these two local spinon singlet states 
are non-degenerate: the energy of the two-spinon state 
differs from the zero-spinon state by the energy cq -|-4u/. 
With non-zero A, one obtains as eigenstates simply two 
different linear combination of these two states on each 
site. The lower energy of the two will be realized in the 
ground state, and the upper energy state has no physical 
significance. Physically, the lower energy state, which is 
neutral (rir = 0) and is a spin singlet, corresponds to 
the local vacuum, i.e. a site with no electron on it. If 
Tir = I, then we must have one spinon on this site, which 
may have either spin orientation. This state has thus the 
quantum numbers of a physical electron. Fixing the total 
charge of the system Q = X^^^r 7^ ^ "^^^^ require some 
number of electrons in the system. At zeroth order these 
are localized, but at 0(tcts/h), the electrons acquire a 
hopping between sites, and one obtains a system clearly 
in a Fermi liquid phase (it is not non-interacting, since 
one has a hard-core constraint in the limit considered). 

The above considerations can be applied for zero or 
non-zero A, and there are no qualitative differences in the 
results (the detailed nature of the vacuum state depends 
smoothly on A, leading to a weak dependence of the 
effective electron hopping on the ratio A/uf) in either 
case. This strongly suggests that A is not a "relevant" 
(in the renormalization group sense) perturbation in the 
Fermi liquid phase. This notion can be confirmed more 
formally by considering the limit of very weak A <^ Uf, in 
which it may be treated perturbatively. At zeroth order 
in this perturbation theory, the vacuum state (on a single 
site) is just the state with zero spinous. Formally, the 
perturbative relevance of A is determined by the behavior 
of the two-point function of the pair-field operator, e.g. 

Ca{t) = (/rT(T)/H(T)/rVO)/r\(0)). (A6) 

Since the pair-field operator creates a site doubly- 
occupied by spinous, the energy of the intermediate 
states encountered in the imaginary time evolution from 
to r is increased by w,y, so that the spinon pair-field 
correlator decays exponentially, Ca(t) ~ e"'*"-^'^. This 
indicates that the spinon pair field is strongly irrelevant 
(formally with infinite scaling dimension). The impor- 
tance of this observation in the context of this paper is 
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that it provides an example in which the spinon pair field 
- which naively has a special significance because it alone 
violates spinon number conservation - is irrelevant. This 
irrelevance is a consequence of strong vorticity fluctua- 
tions, which bind (confine) charge to the spinous to form 
electrons. Since charge is conserved, electron number 
must be conserved in the resulting effective theory. Sim- 
ilar (but not quite so large) vorticity fluctuations in the 
RFL have the effect of rendering the spinon pair fleld 
irrelevant. 



APPENDIX B: ENSLAVING THE Z2 GAUGE 
FIELDS 

To enslave the Zi gauge fields, we employ two sequen- 
tial unitary transformations, C/12 ~ UiV^, with 

/ 00 \ 

r \a:'=0 / r 

/ 00 



The two operators \J\ and \J2 are mutually commuting. 
Roughly, U\ transforms to a gauge in which erf = 1, and 
Ui transforms to a gauge with 'a\ = \. More precisely, 
applying the first unitary transformation, H^i and Hjss 
are invariant, while the fermion Hamiltonian transforms 
to 



where 

2. slave 



a 



1 

2. slave 



(r) 



(B3) 



(B4) 



w-n n -^-rm; 

a;'=0 □(r+w+s'x) a;'=0 

The vortex kinetic terms also transform 

^\'^kin^\ — "-kir. 

with 

\^) = 1, 



+w+x'>i 



— T/ \ — 1 j:, slave/ \ 

a.{x)~^a-CF- (r;n) 



-a;, slave / 



(B5) 



(B6) 



slave J^J ^ -|^^n^^___ 



Simultaneously, the first constraint is rendered trivial 

(B7) 



^l,slave ^ U\CIU^ = n = 1- 
□ (r) 



Further transformation with Ui removes cr^ from Rll^ 
and trivializes the remaining constraint, i.e. 

^2,slave ^ ulc^U ^ = H = (^^) 
□ (r) 



The final transformed Hamiltonian no longer involves any 
dynamical gauge fields (whose state is uniquely specified 
by Eqs. ljB7l IBSp ). and is simply given by 



(B9) 



The electron destruction operator in the Z2 vortex- 
spinon theory, c^a = Wfra with in Eq. 140() . transforms 
upon enslaving in an identical fashion. 



(BIO) 



APPENDIX C: ENSLAVING THE (7(1) GAUGE 
THEORY 



As for the Z2 case, to enslave the gauge fields in the 
U{\) formulation we apply two sequential unitary trans- 
formations, Uab = UaUb, with 



Ua 



- p*Er r' N,V(r-r'){V-S)(r') lY,^ ^, y(r-r')(V-/3)(r') 



(CI) 
(C2) 



where V^F(r — r') = 6rr' and, 
where the lattice Laplacian is 

VV(r) = ^9|/(r-x,) = ^/(r+x,)+/(r-x,)-2/(r). 
j j 

(C3) 

Here we have introduced a "angle" function 0(r — r'), 
to be determined later. The two unitary transformations 
commute with one another. Again, both Hpi and Hn 
are invariant commuting with Ua and Ub, whereas under 
application of Ua the vortex kinetic energy transforms to. 



UaHkinUa — Hkv 



e J — >e J e i 



(C4) 



Here, a* is the transverse part of a,- and a*''""^ satisfies. 



7rn 



slave 



0. (C5) 



Similarly, the spinon and electron hopping Hamiltonians 
transform to, 



U^Hs/eUa — ffs/e| 



ig. i3sl<ii'E(JV) , (C6) 



with /?' the transverse part of Pj and 



V • P"^'''" = 0. (C7) 



Notice that the longitudinal parts of a and /3 have been 
eliminated, and the remaining transverse pieces commute 
with one another and can be treated as c- numbers. Si- 
multaneously, the two U{1) constraints are rendered triv- 
ial. 



-tslave 



UlGfUa = A.^M9."*(r-w) ^ (^g^ 
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gsiave ^ uIQvU^ = Z^r^'^'J^'^'j^''"'^) = 1, (C9) With this definition, one finds 



implying that a* = /3* = 0, and fully eliminating both 
gauge fields from the full transformed Hamiltonian, 



UlHU„ = H\ 



(CIO) 



Further transformation with Uh, which is essentially a 
non-singular gauge transformation of both the vortices 
and spinous, modifies the enslaved gauge fields (so that 
they vanish on the horizontal bonds and are integer mul- 
tiples of vr on the vertical bonds) while leaving the gauge 
fiuxes invariant. Specifically, we require 

9,e(r) = 2^,(r), Vr, (Cll) 
dyQ{Y) = Wy{Y) Vr 7^ a;x -)- w, X > 0. (C12) 

Here 13 j (r) is the enslaved gauge field configuration for a 
vortex located at r = (i.e. determined from Eqs. (|C7p 
with Nr = Sr.o). The transverseness of /3j implies then 



V2e(r + w) = V:.(5y,0"'5j,,_i), 



(C13) 



with an unknown function ipx such that il)x = ^ for x < 0. 
Taking a line sum around the origin requires then 



5ye(a;x + w) = -27r + 2/3j^(xx 4^ w). 



(C14) 



for a; > 0. Eas lC13IC14l are the lattice analog of Laplace's 
equation and the condition that Q jumps by 27r across 
the positive x-axis. These conditions determine ipx and 
hence Q. After some algebra, the solution is expressible 
as a Fourier integral 



e(r) 



e(k)e*'^■^ 



where 

e(k) = - 

with 



27re 



ik-w 



|sin(fc,/2)| 



Asin{kx/2)^sm^{kx/2) + 1 



(C15) 

(C16) 

(C17) 
(C18) 



For large arguments, the asymptotic behavior can be ob- 
tained, 



0(r) arctan[j//a;]. 



(C19) 



for + y'^ 3> 1, with the arctan defined on the interval 
[0, 27r]. Hence 8(r) gives a proper lattice version of the 
continuum angle function. 



(C20) 

.,1 — X, slave/ f\ J z. slave/ j^t\ j n i • 

With a- [r;n-') and a^' (riA*) as defined m 
Eq. ljB6|l and Eq. I)B4|I . respectively. Remarkably, the 
enslaved U{1) Hamiltonian is identical to the enslaved 
Z2 Hamiltonian in Eq. H'^""" = 

We have thereby established the formal equivalence be- 
tween the Z2 and U{\) formulations of the vortex-spinon 
field theory - the unitarily transformed enslaved versions 
of the original Hamiltonians and H in Eqs. (|28I43() 
are identical to one another. Finally, we can verify that 
the enslaved versions of the electron operators in the U{1) 
and Z2 formulations also coincide. From the definition 
of the electron operator in the U{V) formulation, Cro- in 
Eq. H55|) . one can readily show that. 



(C21) 



With the analogous expression for the enslaved Z2 elec- 
tron operator in Eq. HljlO() . and upon comparing the 
defining expressions of the electron operators in the Z2 
and U{1) formulations in Eqs. (|41|l and (|55l) . respectively, 
one thereby establishes the desired formal equivalence: 



APPENDIX D: THE VORTEX-ELECTRON 
FORMULATION 

In this Appendix we briefly discuss a third Hamilto- 
nian formulation of the vortex-fermion field theory. The 
vortex-electron Hamiltonian will be expressed in terms of 
"electron operators" , or more correctly operators which 
create excitations having a non-vanishing overlap with 
the bare electron. In order to transform to this formu- 
lation, we start with the enslaved version of the U{\) 
vortex-spinon Hamiltonian as obtained in Appendix lUl 



H\ 



(Dl) 



where af""^"^ and j]^^"-^'^ satisfy, 

eijdiaf'''"'{Y - w) = Tini- V • d'^""" = 0, (D2) 

e^J^^pf''''^r -w)=t:N,; V • 0'^"'"'' = 0. (D3) 

Now consider the unitary transformation, 

where again V^y(r — r') = Srr'- As is apparent from 
Eqs. (|39I59() this transformation takes one from the 
spinon operator to the electron operator. 



rcr — ^ra- 



(D5) 
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Here ^^(r) = 0^ e^"""* is defined in Eq. ^ and the 
last equality follows from Eq. l(5^ in the enslaved gauge 
with purely transverse gauge field, (3^ = 0. The electrical 
charge density in the vortex-spinon formulation trans- 
forms to include the electron density: 

Ulitijdiaj{v - w)Uei = eijdiaj{r - w) + ncl^Crcr- (D6) 

The full Hamiltonian density within the vortex- 
electron formulation is readily obtained from the enslaved 
vortex-spinon Hamiltonian: 7i„e = U'^iUlTLUaUei- It can 
be compactly expressed as, 

3 " 

+ n^,{aj)+Hr{aj) + Hf, (D7) 

where Ti.y and Tir denote the vortex and roton hopping 
terms, respectively, and are given explicitly as, 

Hv{aj) = -t^ ^ cos(aj), (D8) 
j 

TCrittj) = --y ^ cos{etjdiaj). (D9) 

i 

The fermionic Hamiltonian is 

^/ = -E*-4*.-^ra (DIO) 

j 

i 

The electric field appearing in the spinon hopping term 
yields the same physical effects as the gauge field in the 
[/(I) formulation. Specifically, when the electron hops 
from one site to a neighboring site, the factor e*'^'^ which 
shifts the gauge field aj by tt, effectively hops a com- 
pensating chargon in the opposite direction. In addition, 
this minimal coupling form encodes the requisite minus 
sign when a spinon is hopped around a vortex and vice 
versa. The full Hamiltonian must be supplemented with 
the constraint that W ■ e — N , with integer N. We em- 
phasize that the total spin, 

5 = ^E4.WCr,„ (Dll) 

r 

and electric charge, 

Q ^ J2[^e,,d,a,{r - w) + ct^c^,], (D12) 

z — / 7]- 
r 

are conserved, commuting with Tiye- 

It is of course also possible to pass to a Euclidean path 
integral representation of the partition function associ- 
ated with the above vortex-electron Hamiltonian. Specif- 
ically, the corresponding Euclidean Lagrangian can be 
readily expressed as, 

Cye = iejdoaj + cldoc^ + iao(V • e — N^) + Hve, (D13) 



where the time component of the gauge field ao(r) lives 
on the sites of the dual lattice. In the partition function, 
the vortex number is a continuous field running over 
the real numbers, but the integration only contributes 
when N is integer. To see why, it is instructive to let 
fl/j ^ o-fi^dp,9, and to integrate the vortex phase variable 
9 over the reals. Since the Hamiltonian is 2tt periodic in 
6, upon splitting the integration as dgO = 2-kI + OqO with 
6 = [0, 27r], the summation of exp{i2TT£N) over integer £ 
vanishes unless the vortex number TV is an integer. 

To obtain a more tractable representation of the La- 
grangian, we introduce a Hubbard-Stratonovich field, 
eo(r), to decouple the Coulomb interaction term above. 
Here ieg has the physical meaning of a dynamical elec- 
trostatic potential. In this way the full Euclidean La- 
grangian can be conveniently decomposed as the sum of 
a bosonic charge sector and a fermionic spin and charge 
carrying sector, Cve = Cc + Cf. The full bosonic sector 
is given by, 

= -^[Cj + -^-^(^ol + iejidoaj ~ djao) - -eoeijdiQj 

Zi TT "t^Q TT 

+ i(ao6l - ao)A^ + -Ct, + C (D14) 
with vortex and roton hopping terms, 

Lv = — tt, E cos{dj9 ~ aj), (D15) 
j 

Cr = -Y[cos(A^.y6' - d^ay) -t- (x <-> y)]. (E)16) 

The Lagrangian density in the fermionic sector is, 

Cf ^ clido - ieQ)cr + Hf + ieopQ. (D17) 

In addition to the global symmetries corresponding to 
spin and charge conservation, the full Lagrangian has a 
local gauge symmetry, being invariant under. 

Or — > + ©r, 

a^{r) ^ a^(r) + 9^e„ (D18) 

with Or an arbitrary function of space and imaginary 
time. Because of this gauge invariance we are free to 
choose an appropriate gauge. 

APPENDIX E: POLARIZATION TENSOR 

In this appendix, we calculate the corrections to the 
polarization tensor H^ at 0{t^) including fluctuations 
of both d and 9, for the general case fo < oo. Inte- 
grating out all dynamical fields to 0(t,y), one finds that 
the effective action as a functional of Aj takes the form 

Sf ^S°^ + S^^\ where 

S'l'=-P'^j:[ (a(r,r)C,(r',r')e-^A , 

(El) 
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with the shorthand notation, Ci(r, r) = cos{di9 — ai)rr 
and with the {■)a=o indicating a Gaussian average with 
respect to the RL action. This can be written as 

= -|e^- (^e^^l)--o + e{r?)..o) ^ (e2) 

with 

-^W.„ICy9ik,LJn)], (E3) 

= e'(k-r-^„r) _ g«(k.r'-a.„r')^ ^nd obtained from 
Fx by a; <-> y. Evaluating the expectation value gives 

^ £ L i-'+vUr-r'r]^^ ^""^^ 



+.X ^ y. 

Since we are interested in the polarization tensor, we may 
expand the exponential in Eq. ljE4|l to 0{A^) to obtain 

Hy =nf>,. + n|f , with 

X \ICyGi2+ICxGn\^K,Kj + ix^y). (E5) 

In the limit of interest for the conductivity, |k| ^ at 
fixed frequency, Gu ^ G12, and we obtain Eq. p53|l of 
the main text. 
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